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Abstract 


A data tape is necessary for the storage of a 
systematized collection of physico-mechanical proper- 
ties of rocks. Specific programs permit the obtaining 
of descriptive information on the data - ranges, means, 
and counts. Statistical routines yield histograms, 
scattergrams, and least squares equations. One objec- 
tive is to provide information that can form the basis 
for some degree of uniformity in such research. Ex- 
perience with the programs yielded certain principles 
and changes that would improve efficiency. It is 
recommended that one choose an efficient means of 
data storage, maintain a back-up data source and pre- 
cise records, run all descriptive programs first, 
process as many cases per run as possible, determine 
common scales where necessary for later interpreta- 
tion, and label output meaningfully.. 


l. INTRODUCTION 


Several requests have been made for information on how the 
statistical analyses used in Technical Report No. 2 have been 
handled, especially with reference to the computer techniques 
involved. The information in this supplement to Technical Report 
No. 2 is a brief summary of these techniques. It includes both 
specific examples of computer programs used to obtain the results 
presented in Technical Report No. 2 and commentary on variations 
and possible improvements. These programs are intended as a use~ 
ful guideline for those doing similar research on their own data 
at their own facility. 


2. DATA STORAGE 
SE Ee. Se 


The preliminary work involved transferring the data to mag- 
netic tape for efficient storage. First, the data were recorded 
on the coded sheets described in Appendix A of Technical Report 
No. 2. The data on these sheets was then kevpunched on standard 
80-column IBM cards. A scratch (temporary) tape was then created 
from these cards. The data on this tape was subjected to the 
standard UTILITY SORT routine for the IBM 7094. This procedure 
provided a second scratch tape which contained the information 
sorted according to rock type. The final operation involved the 
use of a COBOL program to create a magnetic tape containing the 
original data, sorted on rock type, in such a manner as was de- 
sired for print-out. Note that the SORT routines often distin- 
guish between a blank and a zero. Hence, it is to one's advantage 
to be consistent in this respect when recording the data on the 
coded sheets. For example, in our case, a three-character field 
was allowed toaccommodate the number identifying different rock 
types. When this number contained only two digits, say 22, it 
was sometimes entered as 022 and sometimes as _22. Either way 
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provide the necessary identification for the corresponding 
data. But since the SORT routine considers these as two dif- 
ferent things, the group of data for rock type #22 on our tape 
is split, an undesirable situation as part appears in one 
section of the tape and part in another. 


Several copies of this tape have since been made because 
constant use eventually results in bad spots that produce READ 
PARITY errors. (One good copy should always be reserved as a 
back-up tape.) Also, it was convenient to have several tapes 
available in order to speed up the computer turn-around time. 
For example, if each of two programs requested either tape 
#544 or tape #630, it was possible to have both programs run- 
ning at once, each using the same information but from different 
tapes mounted on different tape drives. This could not have 
been done if both programs requested tape #544 only. One of 
the programs would have had to await the return of tape #544 
before being run. (Our original tape was created for the IBM 
7094 because that was the computer in use at that time in our 
Computer SciencesCenter. Since then, the CDC 6500 has re- 
placed the 7094 as the principal computer being used. We are 
now in the process of converting our tape to one internal to 
the CDC in hopes that computer usage of the tape will be fe- 
cilitated. Then, because of better turn-around time, it may 
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no longer be necessary to maintain several tapes.) 


3. DESCRIPTIVE PROGRAMS 
3.1 LISTING OF TAPE (pp. 13-15) 


The first step after creating the data tape is to obtain 
a complete listing of the tape. Only in this way will you be 
assured that all the data have been successfully transferred 
to the tape. Such a listing will also provide a useful record 
for future reference. It should be noted agin that one must 
always maintain a duplicate tape as a back-up tape. Accidents 
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or repeated usage can result in a tape that can no longer be 
successfully read by the computer. 


After having obtained a listing of the tape, one may want 
to run programs that will provide for detailed information 
about the data on the tape. Several of the included programs 
are of tnis nature,while others perform certain computations 
on the data. Each program will now be cormsidered individuallv. 
In the Conclusion, I will comment on other similar programs 
that covld be run and on certain approaches or changes that in 
retrospect would have been more efficient and certainly less 
time consuming. 


3.2 COUNT PROGRAMS 


Program + 1 (pp. 16-23) is run to obtain information on 
the number of non-blank entries present for each of the vari- 
ables on the tape. Each variable on each of the eight data 
cards per set2/ is read and compared with a constant, B, that 
is preset to blank. If the variable in question is non-blank, 
i.e., not equal to B, a counter is increased. (Note that the 
counters are initialized to zero as a precaution.) Otherwise, 
the counter maintains its value and the program proceeds to 
the next test, until all variables have been compared. The 
program then reads another record from the tape and proceeds as 
above. This process is repeated until all the data on the tape 
is exhausted. The output from the program consists of a list- 
ing of the variable names and the associated counts or N's. 
(Note that the variables are tested against a blank value and 
not zero because zero is a legitimate entry, i.e., a blank in- 
dicates no data whereas zero is an acceptable value.) 


If such a count is needed only for certain selected vari- 
ables, then one need read only those particular values from 


Vote that the eight cards that are referred to as a "data 
set" in Technical Report No. 2 will be referred to as a 
“record” with respect to the data tape. 
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the tape. (See program ia, pp.24-26). Depending upon the manner 


in which your data is set up on your tape, it may or may not be 
necessary to provide for reading dummy information in order to 
skip to the next record. For example, suppose we wish to con- 
cern ourselves only with a variable - say A20 - on the A card. 
Then in order for our DO LOOP to execute properly, one of two 
alternatives is necessary. We could have the following READ and 
FORMAT statements: 


READ (1,1) A20 
1 FORMAT (10X,A4///////) 
where the seven slashes would cause the computer to skip the next 


seven cards and begin reading again with the first card of the 
next record, or we could read dummy variables as follows: 


READ (1,1) A20 


READ (1,3) B 
READ (1,3) C 
READ (1,3) D 
READ (1,3) E 
READ (1,3) F 
READ (1,3) G 
READ (1,3) H 
1 FORMAT (10X,A4) 


FORMAT (50X,A5) 


where the tape unit we are reading from is Tape 1. The latter 
alternative is selected as safer since the interpretation of the 
slashes is subject to variation depending on the particular 
computer. 


Program #2 (pp. 27-30) is a slight variation of Program 
#1. For certain chosen variables, it is necessary to know how 
many cases there are of non-blank entries for both variables for 
tne same data set. For this program, the required variables for 
each record are read and then each variable for the particular 
pair is compared against the blank. If both are non-blank, the 
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counter is incremented. Otherwise, the program goes on to the 
next test. After all pairs for a data set are processed, the 
variables for the next data set are read and compared. This 
procedure is repeated until all data on the tape is used. The 
output is a listing of the values of N for the various pairs. 
Note that here, as previously, the reading and comparisons are 
performed within the same DO LOOP. Because of the large quan- 
tity of data, it is not feasible or even possible to read the 
data into arrays and then perform the comparisons separately. 


One further variation of these programs may be useful. It 
May be necessary to obtain such information for the variables 


within a particular rock type rather than over all rock types. 
Such can be accomplished by reading in the rock type variable 

along with the data in question and testing its value against 

the vulue identifying the rock type of interest. The sequence 
of instructions would be as follows: 


DO 100 K = 1,2170 


READ (1,1) IRT,A20,Al 
READ (1,3) BB 


am 


READ (1,3) HH 

FORMAT (15X, 13, 10X, A4, A3) 

FORMAT (50X,A5) 

IF (IRT .NE. 223) GO TO 100 

IF {(A20 .NE. B) .AND. (Al .NE. B)} N=N+1 
100 CONTINUE 


where IRT is the variable for rock type and 223 is the value of 
the rock type of interest. 
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3.3 RANGES AND MEAN VALUES 


Program #3 (pp. 31-35) and Program #4 (pp. 36-43) pro- 
vide further information about the data on the tape. Program 
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#3 is designed to print out all values within specified ranges 
for a particular selection of variables and the associated rock 
types. In addition, the program computes for each variable the 
maximum and minimum values, the number of values, and the mean 
over these values. 


At this point, it is important to understand that the data 
on our tape was entered in such a manner that it could be read 
off only in integer or A field format. The variables could not 
be read off as real variables. For all homogenacus computa- 
tions ‘i.e., computations involving only one variable) it was 
not necessary to perform any scaling in order to obtain mean- 
ingful recults. Any scaling could be applica directly to the 
results. So if the entry 432 for variabie B8 really repre- 
sented 4.32 and we got a mean of 397 for B8, we knew immedi- 
ately that the actual mean was 3.97. However, if the computation 
involved more than one variable, it was mandatory that the vari- 
ables be scaled before the computation began. This fact will be 
discussed later in conjunction with several of the computational 
programs. Also discussed in conjunction with the computational 
programs is the reason for specifying ranges for certain of the 
variables. (pp. 7, 9 and 10) 


Program #4 differs from Program #3 in that the print-out is 
limited each time to a particular rock type. This permits the 
comparison, for example, of the ranges of values for a particu- 
lar variable over the whole tape with the range within a spe- 
cific rock type. The actual differences in the program are 
modifications to test for rock type and to reinitialize coun- 
ters. It is, of course, possible that some variables for a 
particular rock type would be blank, and such a situation has 
to be taken into account. It is also possible that the last 
non-blank record read from the tape would be of a different 
rock type than the preceding. In that case, it is necessary to 
provide a way for printing out the last informatior. By insert- 
ing the test on END OF FILE, it is possible to reinitialize all 
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values, go ru: k and perform the tests, and then print out the 
desired inforzation. The test on END OF FILE Guarantees chat 
this last informaticn will not be lost either through exceeding 
the DO LOOP limits oer through trying to zsead past the END-~-OF~ 
FILE on the tape. (Note that the DO LOOP p. vameter goes from l 
to 2170 previously and from 1 to 2171 in this case to allicw for 
this laet possibility.) 


The following are also to ne noted. The manner in which 
the DO LOOP parameter is set up requires that the first time 
through the loop you test IRT(1) against IRT(0). (The IRT array 
represents the values identifying the different rock types.) 

To make this test possible, we equivalence the first element of 
the IRT array, i.e., IRT(1} with the second element of a4 two- 
element array called DUMMY. Hence IRT(1) = DUMMY(2). Then, 
IRT(0) is equivalent to DUMMY(1) which is set. equal to the value 
of the first rock type on the tape. In this way, the first time 
the program reads the rock type, it reads IRT{1) = 3 and when it 
goes to compare IRT({1) with IRT(0),. the computer has a value for 
IRT(O) and it is 3. From there on, all proceeds as usual. 


4. COMPUTATIONAL PROGRAMS 
4.1 LEAST SQUARES 


Program #5 (pp. 24-26) also reads selected variables from 
the tape and compares certain pairs with blank. When both vari- 
ables in the pair for the same data set are non-blank, the val- 
ues are stored in arrays. However, the values in this case must 
be scaled so that the Least Squares Equation will make sense. 
For example, if 429 for B15 represents the real value of 4.29 
and 364 for Al is actually 36.4, this difference in decimals 
must be maintained for the resulting equation to be meaningful. 
In addition, the associated rock types are stored in another 
arvay, aud this information (.i.e., the pairs of points and the 
rock types for which they occur) is part of the print~-out. The 
vaiues in these arrays are the input data for the Least Squares 
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routine. The routine computes the least squares coefficients A 
and B of the equatien Y = AX + B from this data. Note that the 
comment cards in the Least Squares program deck adequately ex- 
plain its source and setup. As meny equations as desired can be 
obtained in a particular run, suhject of course to the storage 
limits and time considerations of a particular computer. It 
should be noted that the N for a particular pair of variables 
must be known in advance. This informaticn is needed not only 
for the dimensioning of the arrays but must also be fed as a 
parameter to the Least Squares routine. Hence, we see the value 
of the Count Programs. 


4.2 HISTOGRAMS AND SCATTERGRAMS 


These Count programs are also necesgary for the histogram 
and scattergram routines which follow. These require the number 
of data points as one of their input parameters. The programs 
are set up in such a way that the first part of both reads the 
selected variables from the tape and makes the comparisons with 
the blank. For those cases where both variables have non-blank 
entries for the same data set, the program stores the values in 
arrays. After the whole tape has been read, these arrays are 
put on the disk and serve then as input for the parts of the pro- 
gram that produce the scatt2rgrams and histograms. These parts 
of the reutines are taken from: "BIOMEDICAL COMPUTER PROGRAMS, 
HEALTH SCIENCES COMPUTING FACTLITY, DEPARTMENT OF PREVENTIVE 
MEDICINE AND PUBLIC HEALTH SCHOOI. OF MEDICINE, UNIVERSITY OF 
CALIFORNIA, LOS ANGELES, JANUARY 1, 1964", W. J. Dixon, Editor. 
Included with the print-out of these programs as they have been 
adapted to run on our computer are the user's write ups. These 
provide information on how +o set up the data cards, and the 
various options available. For example, the program that gives 
the scattergram output could be run just to obtain correlation 
coefficients. 


You will also note that these BMD routines are on a library 


file at our facility and are loaded as part of the user's program 


via the control card 
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LIBCOPY (STATBIN, LGO, BMD2D) 


Hence, the program deck for these two routines is set up as 
follows: 
Control Cards 


7 
By 


Program Deck 


7 
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Data Cards for BMD Routine 


iP | as 
9 
This set up is not immediately obvious from the print-out sheets, 


especially since the BMD routines are not printed out each time. 


Another important point is that these routines accept only 
real variables as input. That is, the data format card is valid 
only in Fw.d format. It does not matter that the variables in 
the first part of the program are written on the disk in A for- 
mat. The BMD routine will automatically convert them as it 
reads them. What is necessary though is that the Fw.d format 
accurately represents the true value of the variables with re- 
spect to each other. For example, suppose two variables are 
written on the tape via the format (A4,A3). Then if the values 
1492 and 765 really represent 14.92 and 76.5, they must be read 
by the format (F4.2,F3.1) and not by (F4.0,F3.0). The former 
will permit the BMD routine to perform meaningful computations, 
whereas the latter will not. 


1492 - 765 = 727 would be an erroneous and 
meaningless computation because the real com- 
putation is 

14.92 ~ 76.5 = -61.58 


The histogram program is also run to obtain a master histo-~ 
gram for each variable on the tape. (For further explanation 
refer to Saction 4.4 in Technical Report No. 2) The program is 
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a simple variation of the one included. As before, one needs to 
know the number of non-blank entries for the variable in ques- 
tion. Then it is a simple matter of reading the tape and trans- 
ferring the values to the disk for input to the histogram. In 
several cases in our computations, upper and lower limits were 
set for the variable and each value was tested not only against 
the blank but also against the desired range. This testing was 
necessary to eliminate bad entries. These latter had the un- 
fortunate effect of altering the scaling on the histogram in 
such a way as to make the graph meaningless. For example, sup- 
pose the range of vhe selected variable was 0 to 20, and thet a 
value of 17.4 had been erroneously entered as 174.0. Then in-~ 
stead of being evenly spread, the histogram would be squeezed 
together on the left, with just one value to the far right, a 
situation which is of little help for interpretation and com- 
parison. In addition, the histogram program can be run to ob- 
tain histograms for variables within a particular rock type. 

All that is required is the insertion of a test on rock type 
similar to the one previously described. 


5. COMMENTS AND RECOMMENDATIONS 


Experience with the preceding prograns has led to sugges- 
tions which we plan to use in the continuing research on these 
rock problems. We expect these suggestions to improve the ef-~ 
ficiency of the operation and to prevent the reader from re- 
peating some of our mistakes. 


RECOMMENDATION NO. 1 

All programs giving descriptive information about the data 
on the tape should be run first. In this way, erroneous or ex~ 
treme entries will be discovered immediately and the tape can 
be corrected or the subsequent programs can be written with tests 
that will screen out such values. This will eliminate the need 
to rerun programs processed before these errors are discovered. 
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RECOMMENDATION NO. 2 

Investigate the possibility of creating a permanent file 
containing one's data. In this way, one eliminates the need to 
mount a tape each time a program is run, and turn-around time 
should be correspondingly faster. We plan to try this approach 
with the next group of programs we will be running. 


RECOMMENDATION NO. 3 

Process more than one histogram or scattergram per run by 
performing comparisons on several pairs of variables at once 
and writing the arrays on the disk in order, one after the 
other. Then include in your deck data cards for each of the 
successive histograms or scattergrams. (Note the repeat speci- 


fications in the write-ups accompanying the BMD routines.) 
This technique has been tested successfully in connection witn 
another problem. 


RECOMMENDATION NO. 4 
Keep a back-up tape. 


j 
! 
RECOMMENDATION NO. 5 

Keep careful and precise records. These may take in parc | 


the form of comment cards in programs. But they may also be 


written notations, etc. Only in this way is one able to go 
back after a period of time and recall accurately what he aia2/ | 
And in addition, only in this way can the information be suc- | 


cessfully handed over to another person. 
RECOMMENDATION NO. 6 
Determine a common scale for each of the base variables | 
in the Comparison Histograms. Use the scale cbtained in the 
Master Histogram for the variable whenever that variable 


2/ 


This is espec.ally important in our situation, because Pur- 
due's computer center operation is in part run for student learn- 
ing in addition to user service. Thus, there is constant change 
in the system and in available routines. This means that 

slight changes may be necessary ay time a program is rerun 

after a time lag. 


appears in a Comparison Histogram. Note that the SELECT card 
we use specifies an impossibly small interval so that the BMD 
program will supply optimum scaling for the vaiues in question. 
This approach was necessary since we had nc advance idea as to 
the ranges of the variables once the comparisons for the pair 
had been made. In retrospect, we realize that we could have 
expanded the program that made the original counts and had it 
give us maxima and minima also. 
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RECOMMENDATION NO. 7 

Use two label cards for the Comparison Histograms. The 
first will identify the variable in the particular histogram. 
The other will indicate the variable of comparison, i.e., the 


variable against whose values the histogram variable was 
compared. 


COMMENTS 

A few additional commerts are directed toward the BMD 
routines. It would be convenient if we could feed these a 
variable N for the number of data points. If this were possible, 
we could eliminate the Count programs and simplify the first part 
of the present scattergram and histogram routines. This task is 
one to be tackled in the future. It may, however, prove more 
advantageous to replace the BMD routines with ones of our own, 
possibly utilizing the Calcomp Plotter for the graphing. 


It is hoped that the above will be helpful to all who are 
involved in similar research. Perhaps some degree of consistency 
and uniformity can be achieved which will make the research un- 
derstandable to all involved and which will permit comparisons. 
It is also hoped that useless and time-consuming duplications of 
effort will thus be avoided. Much can still be done in the area 
of programming techniques. What is included here is by no means 
ie exhaustive or definitive. But it is set forth with the hope 
ve that it is a start toward communication and a help toward fur- 

ther research. 
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2581562 9032112 
25a 
2581012 9032522 
258 
2581 13 028°34 
258 
2581162 903252 
9993 P 1172007 0 
993 
993 
993 
993 
993 
93 
9931112 521 
817 D 29762007 0 
817 
817 
817 
817 
817 
617 
817 
818 D 29762007 Oo 
828 
818 
8138 
818 
818 
8)8 
618 
816 D 29682007634 
816 


616 

8161012 32421 
816 

8161 26 “17 977 
B16 


816 


282 
282 977419 
0205v55 270 
127 20°09 59 26H 
173 
173 62 +26 
0205055 280 
496 
496 80 
0205v55 204 
58 “04 09) 048 
es 
25 028*42 
02 u§ 
070 
02 2 
02 3 264 
02 3 294 
282 
43 *20 994 


1005614 


Sn 49 


04 05 


87 415 


180 


029026 


75 


81 


34 vy 


oY #181005 


77 $19 YIT 


36 


71 
6a 


41 
19 


93 


oe RR eae RGR RRR SC cee RRR 


Program No. 1 
COUNT PROGRAM 


Number of non-blank entries for each variable on a tape 


fee 
nee ee F s eas =, ai oe pieces rc petaerrimete oS af ae ata 


ee 


a ES EIA SAP SSNS SSO DL TEA PET SRE af SPE TE A IE SUIS SN Se 
C + CLEIN & CTEN Cet2 * vO) aI Sive00 
“s T * (OEIN F (UEIN Cetas® EGS tL¥990 
JHA TINUD COl $424000 L © (622N 2 (62d (eed °_ 20) 31 S3¥v000 
_ Do t84N 8 (C82) Cats OltH st 142000. { * (32N = (HIN Ce*s'*® [G41 19¥000 
T #CZ20N ® CLEIN Co® ate eUTH) JT $92000 Tt * (2429N ® (4e)% Ceear® 623041 S$v000 
T @e(92)N & (92)N Ceear® or) al 192000 F [T + (92) = (9E)N (e°S’ 2? wd) aL TS¥000 
LT #¢(S22N. 2 (G49 (otal * wh gd 684000 ; T © (S2)N = (SIN Ce*Ss* LZOISE - S%¥000 
T eCo29N 2 (COL) Cae ae® geo ST 182000 to (929N & (HEIN (ees? * 99) 41 T9009 
T «fE2an 2 (E2Z0N Cet ate 9H) AI S¥2000 LC © (€2aM = (egy Cera * goat. SE ¥000 
T #(222N 2 C220N Cee ae® Gopal 194009 T © (22)N 2 (2e)N Ceter* 9D) 31 l€9000 
LT oCT2Zo™ & (T22N Ceeeai® ob) al S€2000 Too (l2)M = (Telw Cees? *® E3741 S2y9000 
UC #¢(udZin = (OLN Caoea® Eb) 3] Tez000 To > (02) © (02) Feta? 22/351. Té 9000 
CL e(69)N_2 (69)N (273° pep at sss S24000 Toe (6TN 2 (60M Cera ® LD) 41 slvo09 
UL eC8QIN & (A9IN Coane TH) SI {22000 tT > (83DIN = OCBTIN Cee © 6H) ST (l9000_ 
LT (LOIN = CL9IN (eParecto) at s12000 L * (ZDIN ® (LDN Cees oo geal. Sovcoo | 
| ‘T #(992N 8 (Q9IN CoeareLEopar T2000 [ * (9D)N = (90M (eta * Zeal 10v000 | 
| T #(S9IN = (S9PN (eee *UT9I aT 902000 T+ (STIM 2 (STN Ce *3' * 96) 41 Siegc00 | 
! 1 #(99IN = (PSN Coesere 69d aT 194600 : C * C7tim S_fetlIn (eta * Gebal =... . L4€000 
-£ (EOIN = (E9IN Ceeans GOPdt 49000 OO Coe (ETN 2 CETIN Cetst * 99) 3] S9€000 
T #(29)N & (290N Ceeg?e £5) 41 1490090 {> ¢2TIN F (C2UIN Feta? Ew sl 19€000 
T #C(G)M = (L9)N (ees? 997451 ' g§99¢6u00 : _€ * (CON S_CLOIN CeP sa *® 28y at SSE000 
E *Co9m = CO9IN (Hest? SOVST 199000 C+ (OTN = COTIN Cota? Teal Tse000 
tT @#(6S5)N = (65)N Cee ar® #9) 4T _ $9000 T © € GON = ( 679 (ets * S¥ial Svt000 
IT @(8GIN & -(9S)N Cee gr® €9) 47 1990900. LC * € @N_F_Cwlwikts * Rw) aL {seo0u 
LT *(CZGIN S_ (25) (eta? 29pat 99000. tC > ¢ 2M = € LIN Ce? * Lvydal Stcooo 
ww «COT (9G9N @ 199)0N Cees? TO) Al 199000 . ' (+ € GIN = ( wou Cees * OV) al (€C000 
boo) T +(S6)N & (G5)u Cees eta al s€9000 : L* ¢( GIN F_f£ ody Ceca * Gytel_  _.. §=6§2¢000 
L.e(9S)N = (9S)N Caete tl La al TEé9000 —T > ( MIN = ( IN (ets * ve AI (2000 
UL *(ESGIN & CESIN (Coes euts) aT $29000 i [6 € 9% 2 € Cd Cees * Eval slteooo 
T #(2S)M = (25)N Cota? 64) 41 (29060 i : . _€ * € 2in = 1 22N fees 2yyalL CC COTTE O00 
Le CtSyn = (USN Ceeae? Barat. 8t9000 f+ ¢ TON 2 C TN eta? Twat SOEooO 
T +CUSIN & LOGIN Ceeate garal 119000 : (EVOKE SWEET TUSKESHUSE Ui TE KES HOSE YS ytd eU xvid) LyMmMUY ° socooo | 1 
tL *(69IN = (69IN CHa? G4) 41. £09000 ! CL SOUTHS AHS G48 eon tcHtadsentaueie (e%ircGwaw. _ _. £32000 ' 
L.¢(eeIN = (woo (oeas* Gadat LL U109000 : ZJST9* L199 #019 £694 34 8 49499 SGN6 HOS EYE ZAN® in(%*1) Cae 2120090 
Tt *#(Z@ON & CLOdN (eeat *: o3) aT $2$9000 ave xeevusevemvexs® SVOEKZ tose vor eKs Sv evs avteyeos Kul) wry ¥ 2412000 ' 
T *(99)N © 199)N (Hea? Es) al (29000 : SESS OL as OT sO ste d*ee*oa*Saresseatedeia (o*iiGwaa. . _. .. £91000 
T e(S9IN & (GHD Cetas® gaxyat =| §99000 : £3*239°99*S3492ee 4829003 (ESL Ae sy 4€1000 i 
UL oC oem 2 (9IN (aUea** [4 SE 199000 : LSU XS OVE PHY EKG my EXE SMT EVE Soy Px Gtwee xl?) Lymyud c 4£f 000 f 
TL *(EoN & (EOIN (Ceear? ga) al $SS000 Beiezd*odsaswdsegurdgely tr olidvaa . . et~ooo 
T #t2eN 3 _f2e)u (eta? 23)4D 2 _ 189000 63*HI6L9999SDPHISCED*ZIEID (c*t Gvad 590000 
T (TOM 2 CTS)N ftatae? Ga) 4T- $%s000 : CEVOXZ SEV EXZOEV Kol eoyeogeCyizy Eure? xe?) Lvwaud 2 $960000 
T *CO9)N & COMP Cae3r® G3) 4f _T¥s000 i 68*ss*2g* Ie *Setedsew *Zaele (of liQvda 4£0000 
tl #4(oCIn = sans €s000 1 CEVO KS SUE KASS E UE KS HT EK IIE YEXS P eye St HEY? ag th ¥ESTUPKS) LowHug if 4£0000 
UL e(GO2N = CutdN (eeare Es) 4! 1€s000! i Gv atte *SusGysoreertzgely (i sijyaeay 2lo000. 
tT e(ZQ0nN = (2000 Catae® 23)sI. szso00 . 22 oe LWT le* vol oad otoooo 
tT *#(9EIN 8B (G9E)N Cuta® [3a) aI t2g000 “os (win Of $00000 
tT +(CGEIN + {5f)M (me ah? BQ2aT...... 2. ST§000: weevemw Of UC e0acod 
tT *¢ (WEIN F (GE)N (eee * 40)41 ~T1sd00: fs /wl/ou viva Eo00nda 
T * (CEEIN = (CEEIN (eea*® 9G) 4! . $09000 (ALIN NWOTSN awl G €00000 


tT * (ZEN = (2M vega’ oa)st “T0s000; CTIAV EL LALNOZIIGEL® LNGNT AG Id UL ING LNOtinam iI A wvNoIed 


KANE 


090777 
001003 
¢01003 
001007 
691007 
601015 
001015 
£01023 
001023 
00103) 
601031 
001n3T 
002037 
001645 
001045 
001053 
601653 
00106) 
601961 
001067 
061067 
001675 
001075 
001103 
001103 
oOLli? 
001111 
001117 
GO1117 
001125 
001125 
001133 
901123 
00114) 


00114) 

001147 
001147 
001185 
001155 
001163 
091163 
COlIT) 

00117) 

001177 
001177 
001205 
001205 
001213 
001213 
00122) 
0012?) 
N01227 


001227 
001235 
001235 
001243 
091243 
601251 
00125) 
001257 


ane 
yeng 
anh 
Qne 
an3 
ans 
Qgns 
Qn6 
907 
Gnk 


gny 


915 


Q2? 


924 


wR TTL CASON 

FORMATO LHD) 

WOTEE (A, 1840) 

FANRMOT (LM eLOAs]JHVAKLABLE NAME ca De TRA) 

WOLTE EM SGT) nit) 

FORMAT (PHO e buAa LOnhERMEARBIL ITY nee 14) 

wOTTE (64907) NOPS 

FOHMAT HG, . GASP THI Rit SCECTFIC trrvITY, 304.14) 

WRITE C6 .GO 4) WEA) 

For 4d T (pee lu Ae PSHAPPAREN!T SPECIRYIC Gravity. anels) 
WRITE (A 6ONG) m4) 

FORMAL (PHE eT A sy AAMURSRPECIF EEL TyPk OF SPECIE IC UbmivlTye2@xy Te) 
WUT TE CO. IRY NES] 

FORMAT (HO, pA SRHPOSOST Ty espke lad 

WHITE (6,088) N46) 

FORMAT (Pte LAs LORARSURPTIOR eS63074) 

wRTTE (A897) a7) 

FORMAT( HO, 1A e FT MSSLEHOSCOPE HARDAESSe Weta? 

wWGTTE (6.908) n(4) 

FORMAT ( (AN TIAG LTAARRASIVE “ARDO FSS 0G BAe 14) 

WITTE (6,294) WIG) 

FORMAT (pede tude TOMI MPACT TOUOHNE SS eax T4) 

wWPTTR (o,9T5} NCEA) 

FORMAT (pHa e DUA e THE RORPAGATION VEL@LONG wAVES/LAHe 30Ke 14) 

WHY TE (6,931) ®(1)) 

FOMMAT (pH e) CAs TAHPE OR OGATION VELSTRANSVECSE WAVES/L Ane AaXe Y4) 
WRITE (6,927) N112} 

FAR MAT (Hoe DUA ePOHROEFULPIC BASS ING CAPACITY/L OR e 31 ACT 4) 

ah FR CR GPR) N13) 

FORGAT CY HO. DUA ePRHPCTSSUA S RATYS. DYNAMIC) Re AA 16) 

ah TTR Ge, O15) WL G) 

FORMAT (THN sg Lud ePRBEICENDS S MODELS. DY¥RAR ICs: AWMe JPA OTE) 

WOTLE (6,928) NOES) 

FORMAT CYHGePDACSTRMCRULUS UF MEGCTOTEYs SHEAR MODULUS PDYNEMTC/L ABs 
L1LAX—e14) 

WOETE (t&eQ}e} with) 

FORMAT CERO e PUA ePLHSREAM STRENGIH/STATIC ce 39%- 14) 

WRITE (60517) NCI?) 

FORMAT (PHN e DUA ePRHSEE AR STRFNGT® AT SSI/STATTCe tere] 4) 

BITTE (AD) KR) NCTA) : 

FORAAT (1 rHOePuAePAHSRE AX STRENGTH AT SSZ/STATICe 32%01 4) 
eRETE(A.9DG) NOYG) 

FORGAT (HU e 1 UA e T2HPRORPAGATIGN VELSLONG wKUPS/FE EEL so OAK 9 14) 
wHITTE(E,92f) NOPD) 

FORMAT (PHM 9} UA IRAPGCOPAGATION VE LETRANSVERSE WavESSFIFL D9 72x97 a) 
P2ATTE(Re9P1} WOA]) . 
FORMAT YHE PUA sTHOSPECIFIC (AMP TS CAPACITY/FIEL De 29Xe]4) 
wRITE(6,4?7) NA?) 

FORMAT ( {HM e PUA S0HPOTSSON S RAT Tre DYNAMIC/E TELD a Sue 14%) 
WETTEFE,9P3) NOPD). 

FORMAT (LHe) UAG ACHICUVG S MODULUS». NYNAMIC/ELELDe 30% 514) 
WHITE (62974) (7G) 

FORMATO LHe PA e4SHMCPULUS OF KIGIOLTYs SHEAR MADULUS.s DYNAMTC/F Tel 
Del laeTs) 

WETTE (A,49P5) N75) 


5 FOWAAT (LHe TUAsZTRSRE UR STRENGTH/STATIC/E TFL Me a5 014) 


WET TE (H.97K) NIP) 

FON SAT CAM LAs IGHSPFAN STRENGTH AT SSISSTATIC/FLFELD e 26XKeTa) 
WeTTELF,927) NAT) 

FOR AAT (YH LUA s IGHORFAR STRENGTH ut SS2/STATIC/EIFLDs AHKO TSI 
WOTTEIA.9PR) N( PR) 

FARMAT ( {HM e TUAs 36HCOMPRESSIVE STRENGTHe UNCONFINED /LAR oO 24X— 14) 


17 


0012567 
001265 


001255 
COL27F 


003273 
091301 


00130) 
0033n7 
002507 
0021315 
002315 
601323 
961331 
001331 
091331 
001337 
£01337 
001345 


907365 
C0133 


003353 
01361 


00136} 
001347 
002367 
S0i375 
013735 
601403 
G01«03 
00341) 
001411 
00134)7 
0014617 
001425 
001425 
001433 
101433 
001441 
00144) 
001447 
0101447 
001455 
001655 
001463 
0014623 
001471 
00147) 
001477 
001477 
001505 
001505 
001513 
003513 
001521 
00152) 
001577 


WRITE (6,979) A279} 


929 FORMAT (PAG ee LUA CAT HCC MEE SSIVE Slab rwigleg bwrartal AT LCSE/LARLLTOKe 34 


M 
wRETE CAS Fe) A OW) 


? 


O50 FARMAT (pate TCA ee THUCMPRESSIVE STHENGTRO VSPr ural AT LS2/LARe T9Xels 


S3zL FARMAV CHO PGA ee PMCOMPHESSIVE STs hath Totactal 


Q3e 
aw 
(ets 
Q3as 


Ow 


1) 
wQTTECK,UITY NUT? 


i) 

ROTTER SF ,932> N92) 

Fora (phe tude FARMCNULUS 
wOTTS (6.992) 2023) 


t 


AE LSI/LAR,1AXe]6 


OF MER URSATION AT tS) /E He PTA eda) 


Frama t (peep tAeP2AHPCTSSON S RATED AT USIP AE AS he TH) 


ePPTTE (6, Q26} AE 4) 
wUTTE (ASI) AC 3S) 


FOeMaT Cro et A C2BHTCHAG S MODULES AT US) 71 440 Wars Ta) 
FORMAT (PHD. PUASPEROET POR GPRPER STMESS/L ARs Ake 147 


PHITF (OG, 8H) NCW) 


FO ART i pete DA: IHHCCMPAESSIVE STRENGI os UNCONF INET ZF TEL D922 00 Ta) 


KN TS (6,027) N27) 


QT FORMAT ODM. DUA eS IRCCMPHESSIVE STRENGTH TUTAXTAL AT LSUSFIFL Dc) TAs 


Tg) 
wQETE (A, IIE) NOW) 


GUS POM HET OIRO. PUA Se SRCOMPRESSIVE STYENGTHe TayaxTal AT LS@/FIEI Del 70 


Lita; 
AQP TE (A QAG) NE AY) 


QW FOR AET Ve DetuAea INCCMr -ESSTVE STRENGTH TRTAXTAL 41 LSS/PIFEI Oe] TAs 


240 


G4) 


114) 

SPITE C@,9GN) NEAT) 

FORMAT (Cj 4t 9} UAe JSHNCNYULUS 
wPITThE (6.94}) N641) 

FORMAT CimM, 104s 2HHPCTSSON 
wQTEE GA, 907) (4a?) 
FORMAT FHO ep UAaPRHT OUND S 
aRTTE (64449) 0643) 
FORMAT (PA ebuAePAHDSET FOR 
aRTTE (6,944) N44) 

FORALT (LRM) 1 UAs AGHMANULUS 
WITTE (AR .QaR) N44) 

FORMOT CHM) ) LAs 22HPOQTSSUN 
wOTTE(6,.94H) N46) 
FORMAT] Hh eT UACPAHTCUNG S 
WQ2TTE (66967) (47) 
FORMATO LAMe LUA TIHSET POR 
WQTTE (896%) NOGA) 

FORMAT CP HMe 1} VAs PQHMCOULUS 
SHTTR (OR. 849) NEY) 

FORMAT (1HMe ] VAS AZ EP OLSSUN 
WRITE (6,959) NISMO) 
FORMAT (iH TUAe PZHYOQUNG S 
WOTTF (609851) N(5]) 

Fay AAT (UHM eT UAVTIHSET FOP 
WITTE (AWG?) N52) 
FORMATE) HM eg LUA SAGHMCDULUS 
WOTTE (6,954) (99) 

FOR AAT CLHM 9] UAeA2HKHMOTSSON 
WRITE 66,954) NSS) 
FORMAT (THM e TAs A2HYQUING § 
wATTE (A955) (455) 

FORMAT (IHD) UASITIMSET FOR 
WRITE (6,954) (56) 

FORMAT (PHOT OAs ISRMONULUS 


Ab DEFORMATION AT HSPSETELN 2 25AeT4) 
S RATTO BT US) /FFFLD» 32k0T 4) 
MIDULUS AT USTSFTEL Ne VAxe TS) 
UPFER STRESS/EIFL Ne 34K6 14) 

OF NEF GHM4ATION AT S203} AaTS) 
S HaTlO AT US2e%Rx.T4) 
MODULUS AT US? 3axeTS) 

US2 e49K4614) 

OF QCRORMATION AT (!$3e31) A974) 
S RATIO AT US2eWxe le) 
MODULUS AT USZe Wye LS) 
US3e4GAe14) 

OF DFFURMATTON AT USHo3)X0T4) 
S RATIO AT See aRxela) 
MODULES AT US4eIRx 74) 

US4 04GX—T4) 


OF GEFORMATION AY UISO/FIELN ee 2SAe Ta) 
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001527 watTE(&,SS7) N(S7) 


001535 OST FORMAT (HAs LUA se ABHRCTSSON 

€01535 WATTE (6.958) N58) 

001541 B58 FORMAT (He PUA eARHYQUNG S 

00)543 wOTTE (6,959) NI5Y9) 

00155) O59 FORMAT LHD HAST THSET FUR 

001558} WRTTE(A.960) NAD) 

001557 960 FORMAT (LHOe1 Uz s ISHECNULUS 

001557 MQTTE (66941) N(6)) 

001565 G61 FORMAT (LH. 1 UA ARHPOTSSON 

001565 WRITE £64967) N(N2) 

001573 QK2 FORMAT (HN e LUA s PART ONG § 

001573 SOTTE(A,9462) N(453) 

001601 O43 FARA (pHs TAs) 7HSET FOR 

001601 WRTTE G(R ,Q6a) O54) 

001407 QK& FORMAT (CERO g TULA s ISAECOCLUS 

001607 WETTE (6,985) £558) 

001615 QBS FORMAT (LHe T GAs ZARPCTSSON 

001618 wETTh (6.985) (56) 

001693 S46 FORMAT (YAO. ] UAC PAHYOUNG S 

0601623 wQTTE (6,957) NL4T) 

001631 247 Formal (ihHOeludAg I 7THSET FOS 

001633 WRITE (6_,96R) NOAH) 

001637 ORE FORMAT ( {HO PCAC TGHIENSILE 

001637 WETTE (6.989) NERD) 

001645 QAP FARMAT CHO e PUA SGT HACNULYS 
1) 

001645 WRITE (6.979) NK(7N) 

001653 O70 FORMAT ( IRGe1UAs 32HPCTSSON 

001653 WRITE (6,977) NOT)) 

001661 O77] FORMAT Cie PUA S SPRHICUNG S 

00166) WOTTE (629727) NOT2) 

OO [667 Q72 FORMAT (HO | UAes}HMCOULUS 
ih) 

001667 aRJTE (6,972) 0 (73) 

001675 Q73 FORMAT) HN) UA, 32RPOToSOM 

001675 woHTTE (5,974) (74) 

0012703 O74 FORMAT (EHR e Ae IAAMTOAUNG §S 

001703 WRITE (A e975) N75) 

001711 Q7S5 FORMAT (LHe PUA C4THMCNULUS 
1) 

001711 WRITE (6,074) NC 76) 

OO1TIT Q76 FORMAT ( {He} VAs IARPCTOSON 

OO1TLT WRITER O77) WITT) 

001725 Q77 FARMAT( ERD se} UAeJeHTCUWS § 

003725 WOTTE (6 .47H) ATR) 

001733 OTB FORMAT (Hoe) UA LRRMCNULUS 

001733 STF 
END 


001735 


PROGRAM LENGTH TNCLUNT WG Tru HUFFERS 
006345 


UNUSED COMPTLER SPACE 
012500 


S RATIC aT USP/FIFI De 32% 916) 
MODULES &T USO/FIFILM+AZROTS) 
USO/ETEL Dea LS? 


OF NER UeRMATIOR AT USTGF TFL A5A0 Tb} 


S RATIO AF USAETEL Ms 2X0 fa) 
MODULUS AT USZZETEL Ds 324074) 
US3/FIELDe43N 014) 

OF OFF ORMATYOR 


AT WS4/FJELMNe25Ae Ta) 


S RATIO AT YS4/FTELDe 32X07 4) 
MODULUS AT USG/ETFL Os 32x%+T4) 
USG/F TEL 00434%,14) 


STRENGTHeGGXe 14) 


OF ClUFSSELEY NEFORMAHILITY AT USlalOxe le 


S RATIO AT UST (TENSTLE) «28X9T4) 


MOLULUS AT US] (TENSILE) «PRX 914) 


OF CIFNSILE) DEFORMABILITY AT USPe19A +4 


S RATIO AT US2 ETFNSILE) eQRX9T4) 


(TENSILE) sAAxXeT4) 


“COULUS AT US? 


OF (TENSILE) DEFORMABILITY 47 US2.19X%014 


S RATIO AT USD (TENSILE) e2AXeT4) 


MODULUS AT USA (TENSTLE) ¢p2RXs 14) 


OF RURTHURE 942X474) 


19 


eee eS ee a Nem einen: ae 


{ 
L0aD MAD FILE = LUGO 13.27.0l¢ 12/19/69. PAGE 1} 
Fwa LVAN 100 LWA LOAD 13452 FwA LOADER 57303 FwA VYABLES 56604 UNUSED STORAGE 43132 
PROGRAM ANDRESS FILE Common ADDRESS LENGTH ; | ; 
ve 199 SO | 
INPUTC 6445 SYSTEM ; 
SYSTEM "7&6 4"— SYSTEW "=" SCORED 7566 eae ' 
oureTc 1637 SYSTEM i 
FaTaL78 12177 SYSTEw NO ee | 
Stos J a3TS SYSTEM rr ee ; OE Ee Pee ep yn ‘ 
GETAA 13433 SYSTEM 
/BLANK/ 3 0 i 
| i | 
| ENTRY ANOKRESS REFERENCES (RELATIVE) | 
ve 13 | 
TNPUTC 6447 vc 13 1S 17 2) 23 25 27 31 ' | 
—<s5 4 EEE Ste tg SSS GSR? agg 1 35 36 4) °° 763 ry % eel > | : 
53 55 57 6} 63 64 67 71 | 
73 75 77 101 103 195 107 122 ; f 
RET Bp or fees ep ote ag a Gee eee Ty eek ag Oe ‘lie 115 117 12.00 0 129377077 Yas Yer OT 
133 135 136 14) 143 145 147 151 , : 
153 155 157 16) 142 165 . 167 i7! ; 
EE FES ee ae ee ee ge ip Ere gee 173 17s 177 20) ~ 203 205 SBT ORBIT ' ‘ 
213 215 216 22} P73 225 227 231 : 
233 235 237 24) 243 245 - 247 251 & 
ce ee Pee 252 285 257 261 ~~ 283 265 267 271 
273 275 277 30) 303 304 ‘ 
KRAKER 6551 
QBnrRy—— 1S67 a YESS ce pe ; Fe ee ee Spe ni 
SYSTE* 7774 INPUTS 37 727 
OUTPTC 14 1144 
$/VSTEMC "F741 Be Sage ip ee CE ee ge ae ee ny Mga pe 
SYSTEMP 7767 
END 1663 ve 1736 ; 
SToP ih 2 i ees | eee arama Ye Jc Y Van : EES BN EE oe a ae eK We Ie GR Gee or | 
EXIT 7105 1 
AENORAL 7722 INPUTC 40 730 ‘ 
Santee np Sep ee es eee ee LOUTP EEC 15 1148 " ig SR peat pee } 
QUTeTC 1 )A4) vc 1001 100? 1005 1006 1031 1013 1014 1017 ' 
1071 1022 1025 1027 1030 1033 1035 1036 
Fin ee Se ae SE 1641 1043 1064 1047 1951 2052 ~~ ~—10S55"-—-- 1087 
1060 1063 27065 1066 1071 1073 1074 1077 
110) 1102 2205 1107 1210 2213 1115 1126 
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003436 


PRNRRAM COUNT (INPUT, OUTPUT » TAPEL. TAPE ORZNUTC- (IT) 
OLMENSTION N(77) oR (A) eo RAB) 

RATA A/]TH / 

nO 1 Ke) 927 

Niky & 9 

a9 1000 [al,2)70 ; 

GE AN (192) A2M5A199A) 3eAL 7, AL, AB 

FORMAT (36% AG 2K AG AK GATS ON eAZE AK GAC ICN GAD 
PF AN (1e3) B15eB169RAeR25 

FORMAT (PBY AVA 92x, 49444) 

PEAN(194) Cys 

FORMAT (ARK eA) 

REAN (145) 02790239026 U2 409,026,018 

FORMAT (ALK eA eo TAX 9 Abe 3X Ae OX Ab AS AG OAK, AP) 
PEAN (10H) (RIT) 913198) 

FORMAT (BA10) 

PEAN (197) F3ReF LIF AT FSF IL FSoF 12 

FORMAT (YON AS 5 B3yAG eo) 1X9 Abr AF, 15K, AGA) 
PEAN(196) (WP(T) e123] 08) 

REAN( 168) H2leH? 

FORMAT (21X 9A4 948K 9A3) 


TF AL NEB) N( 1) 2 NO J) 6 2 
TF ¢ 0?,NE.B) N( 2) 2 NC A) © 3 
TF ¢ F3,NE,A) N( 3) 2 N( 3) 6 2 
TF ( F4.NF eB) N( 4) 2 NOC 4) ¢ 1 
IF( FS.NE,A) NE 5) & NO 5) @ 1 
IF( ABeNF.R) N( 6) 2 N( 6) © } 
TF ( H7,NE.B) N( 7) & NC 7) @ 1 
IF( BANE.) N( 8) = N(R) @ 3} 
IF( 09.NE.@) N( 9) 2 NI 9) + 1 
IF(F10,NE,8) N(1O) = N(10) « 1 
TFCFLY NEB) N(TR) = NGAL) ¢ Y 
JTF (F1A.NE.f) N(12) 2x Ni12) + 1 
TF (A13,NE,F) N(13) 2 N(I3) » } 
IF (C14.NE.8) N(14) 2 N(14) © 1 
IF (RIS eNE.B) N(15) 2 (15) o 1 
TF (R16.NE.6) NO16) = N(16) © } 
TF(O17.NE0R) N(17) = N(IT) © 1 
TF (NDANE.B) N(18) = NCLB) + J 
TF (A}9,NE,B) N(19) = N(19) o 1 
TF (A20eNFE.R) N(20) = N(20) ¢ } 
IF (62) NEB) N(21) = N(2}) © } 
VF (N22? NFA) Ni2Z2) = N(22) + } 
TE (N24,NE,B) N(23) = N(23) » } 
TF (N2&eNEoB) N(246) 2 N(24) + } 
TF (R25eNE 2B) N(25) = N(25) + 1 
TF (NAALNE A) N(26) |= N(€6) 6 1 
TF(FO7.NE.B) N(27) = N(27) @ 2 


CONTINUE 

watrE (4,909) 

FORMAT (}H}) 

WRITE (60999) ((KeN(K)) 9K) 927) 
FORMAT(TH 92HN(,1294H) & 914) 
stop 
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' 14626021.C™ 3.566 MwDeSFC. 
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variables in the same data set for 
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000637 WRITE (65.5) SSIUCT) eo TYPEL(T) sB1SAL(T) oA1H1S(1) olzleJl) | 
000660 CALL LSTSQ(J] sNORDReB15a} eALBIS,COEFF) 
000KK4 WRITE (456) IDEN 
000672 WRITE (Hot) (COEFF (T) oT 81 eNROR) 
000705 READ(5.270) IDEN 
000733 WRITE (8.6) J2 \ 
0007) WRITE(605) (S2¢T) oTYPE2Z(]) sBISAL7T (13 sAL7TB15(1) oI) eu?) 
000742 CALL LSTSY(J2sNORDR-B19A1 79 Al 7BISSCOEFF ) 
0007446 WRITE (4.6) IDEN 
000784 WRITE (Se.7) (COEFF (T) 9 T@1eNROR) 
000747 READ (5e20) IDEN 
000775 WRITE(6.4) J3 
001Nn3 WRITE (605) (S3CT) eTYPF 311) eBISAI9(I) ¢A19B15(1) ol#19J9) | 
001096 CALL LSTSQ(J32eNORDReBI SAI 9S ALIBISS COEFF) 
001020 WRITE (4.6) TOEN 
COLON WRITE (4667) (COEFF I(T) os TS1LeNRDR) 
Oo01NSt READ(5—92?0) INEN 
001067 WHITE (6.4) J4 
001065 WRTTE (605) (S4(T) oTYPE4(T) #BISAZ0(T) eAZORIS (1) 9 12) 9 UG) 
001106 CALL LSTSO(J4eNORDReB1SAAN  A2OBISe COEFF) 
091112 WRITE (604) [NEN 
001190 WRITE (607) (COEFF (1) oe T=L9NROR) 
001733 READ(5e20) INEN 
00114) WRITFE(A.4) JS | 
001147 WRITE (665) (SSCT) e TYPES (1) 0B15D02(1) 9926)5(1) oI#19JU5) 
001170 CALL LSTSQ(J5eNORDRsB19U2¢D2H159CUEFF) 
001176 WRITE (406) IDEN 
0012 — WRITE (H07) (COEFF (T) eT 21 sNROR) | 
001715 RFAN(5,20) IDEN 
001223 WRITE (604) J6 | 
001791 WRITE (605) (S6(Ti e TYPFO(T) #B1509(T) eNG815(1) o lel 9J6) { 
001252 CALL LSTSQ(J6eNORDR2B15099098159 COEFF) 
001956 WRITE (AH) INEN 
NGIPKS WRITE (He7) (COEFF (YT) oe T@1sNRODR) 
001277 REAN(5920) IDEN 
0013205 WRITE (606) JT ; 
001313 WRITE (665) iS7(T) sTYPET(T) sAI5F10(T) sFINRIS(T) oL=leJ7) i 
0013234 CALL LSTSQ(U7TeNOROR, BIS" 10¢F L0B1S9COEFF) } 
001240 WRITE (606) IDEN | 
001346 WRITE (607) (COFFF (1) eo lT21leNROR) | 
001941 REA (5920) INEN 
0012467 WRITE (6.4) JA | 
001375 WRITE (605) (SA(T) eTYPES(T) sBISFO(T) oF 4B15(1) ol el eJA) 
0014146 CALL LSTSQ (JA eNORDReB1S5SF49F 48159 COEFF) | 
101422 WRITE (6.6) IDEN i 
001420 WRITE (697) (COFFF(Y) © f=1leNRDR) ' 
001443 REAN(5e20) IDEN 
00145) WRITE (606) JO 
001467 WRITE (605) (SOT) eo TYPEI(T) sBISFLI(T) oFLLBIS(I) el #109) 
001500 CALL LSTSQ(J9eNORDReBLS" LL oF LIBISeCNEFF) 
0N15N46 WRITE (666) IDEN 
001512 WRITE (H.7) (COFFF (J) se IT3]) eNROR) 
0014896 REAI(S,290) IDEN 
001523 WRITE (604) J10 
00184) WRITE (405) (SIC(T) os TYPELOCT) sA1SH2Z)] (1) eH21815 (1) eleledlon 
NOLSA? CALL LSTSACULOsNORNRGALSH2Z1 H2IRISeCOFFF) 
001544 WRITE (406) IDEN 
N01S76 WATTE (407) (COEFF (CY) eI] eNRDR) 
001697 20 FORMAT(1RA4) 
O01497 4 FORMAT (LY H1e///10Xe39HSTRAIGHE LINE LEAST SQUARES FIT THROUGH» 14e1* 


lsL2HPOINTS (XeY) ///) 
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whee reece eae cere ec eee Peeereresscccceteseseeccsecasesercenetoeseel ct S050 


001657 § FORMAT (A(SXoAl eo LHe T30 1H eFSefe LHe eFS.051H))) 

001607 6 FORMAT(1HOs18A4) ' 

001607 7 FORMAT (1HOSSX eo 3HY wef Rede3H * oF Re 3e4h @ X) 

001407 sToP 

90161} END 

PROGRAM LENGTH INCLUNING I/0 RUFFFRS 

020627 

UNUSFN COMPILER SPACE 

022300 : : 

SURROUTINE LSTSQ(NIMeNORDROX oY eRETRN) 

000010 INTEGER RyRR LSTS0002 

090010 REAL MATRX(20920) #YSUM(20) sRETRN(20) LSTS9003 

000010 OIMENSTON X(500) 9 ¥ (500) 9SX (40) 
c LSTS§0005 
Cc mebee tence alee eea sec tet era Oe sateen ssreesee rere esesssosssesscoeeL ST SQ006 
c LSTS0007 
c SUBROUTINE LSTSQA LSTSO00R 
c OFCKS. YVSED LSTS0009 
c LSTSQ LSTSAQ919 
c LSTSN011 
c PURPOSE LSTSQ012 
Cc LSTSQ COMPUTFS AN N THeORDER POLYNOMIAL LEAST SQUARES LSTS9013 
Cc FIT THROUGH m POINTS. LSTSAOL4 
Cc LSTSN015 
c USAGF . LSTSN016 
c CALL LSTSQ(NUMsNORDROXe YeRETRN) LSTSN017 
C LSTSQ018 
Cc DESCRIBTION CF PARAMETERS LSTSN019 
c NUM ~- NUMBER OF POINTS USEN. LSTSA020 
c NOROR = ORDER OF THE DESIRED POLYNOMIAL FIT. S7S$2021 
Cc x = NUM=VECTOR OF X=VALUFS. LSTSN022 
Cc Y ~ NUMeVECTUR OF CORRESPONDING Y=VALUES. LSTSQ023 
C QFTAIN = (NORNRel) -VECTOR OF COEFFICIENTS OF NFSIRED LSTSQ026 
c POLYNOMIAL ORDERED FRom COEFFICIENT OF CONSTANT LSTSON25 
c TERM TO COEFFICIENT OF HIGHEST POWER TERM, LSTSQ026 
C LSTSQ027 
Cc REMARKS LSTSQ02A 
Cc THE CURRENT NYMENSTON STATEMENT IS SUFFICIENT FOR mORAR .LELSTSQO29 
¢C lg AND NUM eLEe 2006 IF THE USER wWISMES TO CHANGE THE LSTS0030 
c NTMENSION STATEMENT TO SUIT HIS NEEDS» THE FOLLOWING FORM LSYSON31 
c SHOULD BE USFN-- LSTSN032 
c LSTSQ033 
Cc TYPE STATEMENT LSTS9034 
c REAL MATRX (NORDR*+1] sNORDR+1) s YSUM(NORDR®)) sRETRN(NORDRe]) LSTSO035 
c OTMENSTION STATEMENT LSTSN036 
¢ DIMENSION xX (NUM) 5 ¥ (NUM) 9 SX (2@NORDR) LSTSQ037 
CG LSTSNO3AR 
¢ Me THON LSTSQ039 
C & STANDARD LINEAR LEAST@SQUARES FIT TS APPROXIMATEN Ay LS1S§0040 
C CALCULATING THE COEFFICIENTS OF THE NORMAL EQUATIONS aNn LSTSQ04) 
¢ THEN SOLVING THE NORMAL EQUATIONS USING TrE GAUSS= LSTSQ042 
Cc ELIMINATION PROCESSe ORTHOGONAL COEFFICIENTS ARE NoT USED, LSTS2043 
C FOR REFERENCE<~= LSTSQ0%6 
Cc MCCRACKENe 0,Ne AND DORNe WeSee #NUMERICAL METHODS ANN LSTSQ0465 
C FORTRAN PROGRAMMING# JOHN WILEY * SONS(NEW YORK==1064)9 LSTSQ046 
c PP, 231-2439 PP.262-275. LSTS0067 
C PROGRAMMER=<RICHARD F. PUK LST§N048 
Cc LS1S0Q069 
c 
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000010 


a00012 
000013 
000018 
000016 
000Nn27 
000010 
000091 


000035 
000037 
060040 
000082 
n00nKs, 
0000565 
00007) 
600072 


000076 
000077 
O001n0 
0001n2 
OO01NS 
000112 
000112 
000117 


0001946 


000195 


000136 
000120 


000132 


000133 
0001433 
0001356 
000147 
000145 
0007140 
C00142 
000145 
GG0147 
000171 
000175 
000203 
00077) 
000213 
000216 
000220 
0002230 
000222 
000240 
000247 
00074} 
000243 
0002c0 
000265 
000257 
000771 
000277 
0004n} 
000In2 
0003s 
000306 


WROR ® NORDRe] 
C 7ERO MATRICES 
DO 16 Tw) sNROR 
YSUM (LT) ee 
NO 10 JeleNRDR 
In MATRX(I,5) 20, 
RR = 2eNQROR 
D0 11 FarsReR 
11 SX(r) a0, 
C CALCULATE Sims 
NO 20 wep eAR 
09 20 ImjsNuUM 
20 SX(K) m§x(K) OX (1) OOK 
NO 25 Kay »sNORDR 
00 2R Tet .UM 
25 YSUM(Ke}) & YSUMEKe 1) (XL) aK) OY (Ty 
NO 2619),NUM 
26 YSUM(1L) = YSUM(}) ey(7) 
C INSERT IN MATRTA (CHANGE FOR DIFFERENT NoRMaL EQUATIONS) 
DO 30 Kat wROR 
DO 30 Ia) eNRDR 
IF (1ak~2)27928¢27 
27 INSERT & PeKed 
MATRX(Key) = SX(INSRT) 
GO Ta 3n 
28 MATRX(K.~) 2 NUM 
39 CONTINUF 
€ SOLVE NORMAl EQUATIONS (MATRIX) 
Cc G4&USS FILIMINATION SFQUENCE 
LITEs2 
09110 Ja) sNRpR 
JI ede] 
TF (JI=NRDR) 60605 
LITEs) 
CONTINUE 
LaJ 
00 12 TeJdleNROR ' 
GO TO (12917) oLITE 
17 TF CARS (MATRA (19.1) ABS (MATRX(L9J))) 12612911) 
YW) Lal 
12 CONTYNUF 
TF (L359 13914973 
13°90 15 KeJsNRDR 
SWITaMATRX(JsK) 
MATRX (JyK) SMATRK (Lyk) 
16 MATRX (1K) aSWIT 
SWITeYSum (J) 
YSUM (J) aYSUM(L) 
YSUM(L) wSwIT 
16 CONTINUE 
09110 TaJteNROR 
GO 10(110923) oLITE 
23 EMBMATRE(T95) /MATAX (ed) 
MATRK (IT. J) BO. 
10 18 Ke) eNRDR 
CBEMOMATRX (J 9K) 
MATRX (TK) @MATRIK( TK) oC 
15 CONT YTAUF 
110 RET ye AeM CL) EM@ySUM Cy) 
RE TRN CNR )#YSUM(NRNR) /MATRX (NROR eNRDR) 
INBNRND 6} 
NO 4) NIX#2snROR 
MIXSY NaN TX 
MINSMI Xe} 
SUMAsD. 
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BMD PROGRAM TO PRODUCE HISTOGRAM 


ee 


SUBPROGRA LENJTH 
090052 


TPwb2 i 


FUNCTION ASSIGYMENTS 


STATEMENT ASSIGNMENTS 


10 = 900057 1é = ovonty 1% - ovonle 1b = 00016 ‘ 
16 = 000014 16 - o00n2e \y = yvonne 22 - Ou0n2¢ 
24 = 1000e6 28 = 000n27 «0 = 000030 49 - 000043 


BLOCK NAMES 4ND LENGTHS 
VARIABLE ASSISNMENTS 


¥ START OF CONSTANTS 


{ 
000040 ! 
1 
STAKT OF TEMPO4ARLES | 
000050 | 
START OF TNITAECTS 
000052 | 
UNUSED CUMPTLER SPACE | 
022600 
PROGRAM BMNSU CINPUTS 1 eCUTPUTS] e FAPESRINFUT, FAPE ORI TPUTs Tapa) ame y0N0N} 
Cowon GENERAL PLOT wlTR H1S1UGHAM AjGUEl IR, 1964 amMOgnn0una 
¢c HEAL Tm SCIENCES GUMPUTING FACIL TIN, “ILA awmhagn0002 é 
c 6400 CUNVENSTON BETTY BFA SON aMOSIHVONSG 
000003 UIMENSTON FG (999) eX (15000) #SYM( ED) ONA (TS) eHE ADT 7) gv4UT(SU0) eXMA(KQ = AMNRDADONS 
10) o Li 15) onXA(LG) AVE SL 917) AMDKDOO0NA 
c aMDsrn0UN7 
govuud3 COMMON XoHF AleX TeSYM oe amGaonduna 
060003 FAVIVALENCE (FGeAMT) ands n00R9 
000003 TYPE ENTEGEN AF232A1 239 C123 L239 I YDE SKY QmMOanQodin 
C aMDSINV0TT 
00v003 110 FORMAT 7OHTSMDNSD GENERAL PLOT @ INCLUDING 4ESTRugaw © VERSTONW UB  AMDSDNUDTA 
1 AUG, lee 1964 / amdeagn0013 
1e04 HEALTH SCTENCES COMPUTING FovleilTY,UCLAas/) RMDSNHVU LS 
Cc amDSNVOLS y 
000003 AL23a (¢onF INISH) aMOspn0014 nai 
090005 nleéd=(*6HPRUGLM) AMON N0U17 ; 
000006 Cl238 (*HHCRSVAR) aMDKINOUIA 
000010 D1l2@3a(*6NSELtECT) am05500019 
000011 i) RE AD ( 5011) TODE se SAME VY PNP ENG ONAN JONTYAN, MT APE PNCARD amMD5i)000¢ 
000037 206 FORMAT (OSHOCUNTROL CAKNS INCORRELTLY ARUERED JQ PUNAUFU) amDsonud?) 
000037 IP (TOVEALES) 2€00e2u1 6200 amdapnvdep 
000062 202 wRITF ( Fe dU) amdann00s3 
000046 20) Stor RMOSDNVUPG 
000050 2u0 IF (TOVESH1E3) 20252055702 QMOK0N00AG 
000052 203 NTAPE RMT APE BMDGINVUA4 
O00uuss IF (NTAPE cE 200) Ni APES aMOG NUOAT 
000056 ING FE INVM(NVO5UL)) S09, 20C e202 AMG DNVOAR 
000062 309 TE CINP OL) ENP 5007) 120862020202 amMdan0Ury | 
060070 205 IFC CWVENADN) @NPH1SNNY) 2V6s200920E am0snn00 39 


000075 206 CALL VFCHCKINCARL) aMD08700031 


Si a aha a rar hr ae 


000077 
009103 
0600125 
000130 
000136 
000140 
e001S2 
000154 
900156 
000162 


00U166 
000177 
000201 
000206 
cdo2le2 
000214 
000227 
a00ze1 
000223 
000226 
000226 
000227 
000231 
000233 
000235 
000237 
000266 
000253 
000255 
000200 
000261 
000300 
000302 
000306 
000310 
000322 
0003z6 
000330 
000332 
000333 
000335 
000337 
000340 

000357 
00036) 


060364 


000370 
000372 
000373 
000376 
000400 
000402 
000403 
000404 
000406 
000410 


aaan 


co) 


211 


9948 
999 


Crs 


ORI TE ( bello) : 
WRITE ( Be C19) SAME gNV, the NG ENAULg VENA SNC ARD 
NTOT eN@ ON yen ; 


LF UCNVONADD) @(NVONATDRSA]  ) LeQurscud 
MCARD ANC ARN 

RE au ( 
NCAKOSNIOTs ve 

DO 2ti LeleVCARD 
n(L) 206u 

NO 3 [#lenp 
Hehan READ TN THE RA DALLA Amd 'KONGPUSFE [AE 4aTorX 
REAL ( NTAPE SF) (Amalia) e vel env) 

OO 3 Galenv 

KANPe JoNPOT 

AUK) BKYALY) 

IF (NTRAN) 20292792) 

CALL THANS (th oNVeNTi2dis) 

LFINV) 202,6€02,22 

NP VSNP 

IF (NAQD) 9949994995996 

NV = NV ¢ NADD 

| 

DO 63 [BleNv 

KMI (1) 899999999, 

XMA11) Se99QI99OQY, 

VO 66 JBKeNPV 

XMI CL) 8 AMINL IX (J) oxMI ery) 

AMA(T)& AMAAL(KX( J} 9 iMAC T)) 

KEK+NP 

NP VSNPVONF 

vO 5) JJal Nb 

ME AD ( 59104) TADE SINT ENL ONL NY 9 Fv 

IF (TQUVE*0123) 202,2v9,202 

1F (NH® (NH#6) ) 21 99eN2 9202 

NHZNHaB 

RE au { 

IF (NC) 202_209R 

NNC®(NCOH) / 7 

IF (NNC #2) 9e%9202 

NG2s0 

VO 150 [ale vnc 

NGLENG24} 

NG@sNGerT 

RE AU ( 59195) TODE® (NACA) e SYA) pUan 3] gp Ntur) 
IF (TODE=C123) 2021 5UerNe 

CONT INUE 


59302) (FG(L) cTz} eteadny 


50102) (MEAD (L) o PEL ot) 


BATTELLE COMMENT 
FILL CUMPUTER WORD wlth INPUT SY MeL 


TF (NC. SE.) CALL FIT LESYM, NU? 


XMAABHA¥IISOINYVYY, 
XMENRGIY9IOIYD, 

LE (NC@1) 205) 1912 

JaNnA()) 

XMAKRXMA CJ) 

AMINSX 4] (0) 

a0 TO 14a 

00 13 Ial>6NC 

JeNA(T) 

XMAAS AMAXL(XMAKSAMA LY) ) 


61 


RMN5HNU032P 
amdsonu033 
5405970026 
amM0cDN00 35 
aAM0aNN0026 
4m08900037 
am)enn003a 
am0cnnd0%0 
AMOS N0060 
8M0G)NU06) 
amd«w n0062 
AMDNGINVDS 2 
RMOA9N0 064 
aMDRNNVNES 
AMDRONNCEA 
QMDENNUORT 
QmM0s0nU06EA 
8M0S0N00469 
RMDSHOU0SO 
RMDS 100051 
aMDRONU0S? 
8M0R000053 
aM05NN005e6 
amOspn0085 
amMOgnn00SA 
aMDC000057 
aMDconuUSR 
RMOENNUDRA 
aMNS0N0N60 
aMDSNAVOAL 
QMDRNO0UH? 
RMDENNDOKR 
RMDENHVUKS 
RMNS5NN0045S 
AMDSNNUURS 
AMDSNHOUDA? 
RMDENADUKS 
RMDKINDOKY 
aMyann0o7n 
aMOGKNNU0TT 
AMOKNAVOT?A 
AMDSONVOT? 
BMORNHNUOTS 
aM0K0N0076 
amMds9n00746 
AMOGHAUVUTT 
amMO5990078H 
aM0e900079 
amoang0gan 
AMOgNNNUS] 
amdaynvugs 
AMOR NNUUNR 
AMDEINVORG 
aMDENNVOAaK 
AMOR INUOAK 
RMDKNNUYUR? 
ampe 00088 
amdan00a9d 
9m9490009n 
aM0aynV09) 
mDENNU0G? 
4m0a0N0094 


Sita att Sei 


rete Sena C cemeanaeneneninanicnasiimemmsiniammensaminil 
ee ee en eG Oe ROR NO Og A Ae AOTC FP Ae RPT eRe ae TT een eA ante mae —eaerarercemnte 


i 
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0004)34 
000623 
000428 
099426 
000434 
000436 
000440 
000046 
0094647 
000463 
000467 
000502 
000503 
000520 
000524 
080526 
060530 
000532 
000533 
000535 
000541 


00056) 
0005646 
000577 
00060] 
000603 
000605 
000606 
000610 
000612 
000613 
000615 
000621 
000633 
000634 
000635 
000636 
000637 
000646 
000647 
000652 
000056 
000655 
000661 
000665 
000667 
000675 
000677 
000701 
000702 
0007046 
000706 
000721 


000726 
000727 
000742 
000TSs 
000746 
000753 
000756 
000764 


13 
i6 
lo 
05 
26 


242 


25 
26 


23 


15 
16 


31 


32 
33 


20 


268 


293 
2846 
2a5 
Cc 29 


29 


35 


RMIN@ AMINT(AMTNG AMT (J) ) 
NP Y¥SO 
00 65 Jel]eNc 


NXA(T) SNK( 1) ONBDeNP 

NY YENY @nNPeNP 

LF (NL) 2323924 

wRALE ( 6ellN) 

1F (9=NC) 296292492249 

WRITE ( He lOBINYS (NA(T) 928359) 
WRITE ( &.11?) 

WRITE ( Ge LIV CNX CTP ed Bho NC) 
GO 10 250 

wRITE¢ Be lLVAINYS(NA(T) 9 bel ent) 
WRIIE ( #2112) 

NO 26 IalsnP 

MYZNYY?] 

yYux iy} 


00 25 Jalent 
MXBNAX (J) OF 


25) eX 14x) 

WRITE ( Feb U6) Vo (Z(K} NBA g NC) 
WRITE UUTPUT TAPE 62310 

WRIEIE (697777) 

WRITE ( 62403) (HEAD (I) 9181 eNH) 
NNO SEN 


YMAAsX4a (NY) 

YMINweX4] (NY) 

nO lo TeleNr 

MYAnYYel 

ysx (my) 

pO 15 J@leNnc 

MKBNXX (J) oT 

Z( J) eX OMX) 

CALL PLUTH (Ye YMINY YMAXK go XMINg KMAK ONG  NNP) 
IF (NNP) 31932932 

NC#e] 

GN TN 33 

NC #U 

CALL PLOTRIVeYMINe YMaAe AMI Ny XMAK SING »NNP) 
60 TO 30 
NYTENY*#NP 
NYYSNYToNPol 
IF ONES 29929028 
WRITE ( 
NNCB(NPO9)/)0 
NG2@BNYY1] 
WRITES fe 
DO 285 {weleNNC 
NGISNG2¢] 

NGCENGEO1G 

1F (NYT#NG2) 283,244,284 
NG 2snyf 
WRITE ( 
CONTINUE 
WRITE OUTPUT TAPE &.!210 
WREVE (osrv777) 
wRITE( 

XMAABK MANY) 
XMENSKME CNY) 2000005 

TF ( (XMAKOERMIN) JF N@= 360) 369349 35 
FNa@(XMAX@XMIN) £300 
WRITE ( 

g0 TO 3 


ello) 


1Oo7) Nv 


he 1€5) (XC) ode Volo NG2) 


Be lU3) (MEAD(T) 95] 9NH) 


6e LU9)FN 


62 


(OXTOTNG ) 


(ORTGINA ) 


gmDROHUOC: 
AMDaGNNO0SS 
AMOR ONVO9K 
RM08H00097 
BMDRI0U09R 
amos9n0099 
amdenognicon 
aeDdcdnvini 
AMDSONVENA 
940GDn0103 
aMORNUI Ne 
amM0enn0lns 
amds9nb196 
am0snn0lo7 
ampbcOnuina 
andannvlog 
aM08200115 
am0snn011) 
amM0s00112 
BMOCIOOL1A 
aMOSOnbI Le 
aMOannvlis 
aMDRDNUIIG 
BMOnNAOLIT 
AMDaNNULIA 
amM0s5on0819 
ReDRNHUIAN 
amMDann0ley 
avDann0l2e 
ampsonula3 
amDspn0lee 
aMDSD00125 
RMDGINVIPA 
amM06190127 
aMOsonulas 
am0s0n0194 
gmMDc Ondo 
am0anndl3) 
amocta0l3> 
3M05900193 
AMOR IN01 4 
AMDNSIOULIS 
QMDE1I)N01 36 
aM0n9NU137 
BMD50NU1 9A 
avo 7n012y 
AMOENNULED 
amDsnn0lal 
PMoaNNULA? 
aM059001463 
BMOSDNUILOS 
aMosnnvles 
aMOaNnO lon 
AMDS9N0167 
aMOgNNOL 4H 
aM05000169 
amMDannulsa 
8M08900151 
aMOGHNULS> 
BMOSONUILAR 
RM08900154 
QM0S900155 
aMDesOULSS 


000765 “3 CALL HISTINYYsNYTeUM TN a XMAA GFN grih) aM0g0nU167? 
00077} 50 CONTINUE BMD8)00188 
000776 GO To 5 awOsn0ulse 

¢ i 5 

000774 TO] FORMAT (2AGeI3e pool aeleeadXeTae2ied) a oeeeriay 
000774 102 FORMAT (8A10) AMDR DOULA? 
00CT74 103 FORMAT (JHOs20x%_HA)N) audson0ias | 
000774 1717 FORMAT (LHD) BMDSNNDI AS 
O0077s 10& FORMAT (AGe2t leo t2eI2.F 11,0) AMIGNNUILGS 
c HATTELLE COMMENT AMOS 0156 | 
i c FORMAT 10S CHANGEU FRUM (7 label 3) 946) BMDRINUIA? | 
c IN ORDER TO COMPLIMFNT oOvEt BATTELLE SUmnROoUFiNeE Ca; Leo FILL. amMOannvles i 

000774 10% FOAMAT(AG! F(19s Rye OX)) BMDENOUIAG 
| OuU0776 106 FORMAT(IM LO(FLUe4.)X)) aMdsnn0170 
000774 107 FORMAT(LH 23H HISTOnRAM OF VARIABLE 13//) AMDSN0171 
000774 108 FORMAT(14H HASE VAPTAGLE >» 3xe16H URUSS VARTASLEC/OX, 1O([398X)) amM0cn0) Fe | 

| 000774 109 FORMAT(LA 954m THE VaLUP GAVEN FON THE INTEavaL atnty IS TOO SMALL aMDRNNOLTA 

| Ll, 733K RA ONEW VALUE SF ll me 2eng HAS HEEN SUDSTTIUTFUYs sy 9mMOn9n0174 
000774 LLL FORMAT OC SX 9S (136K) ) aMoapndl7s ; 
: 000774 Ll2 FORMATCIA ) ‘ ' 
| 000774 125) FORMATCLH 10F 1104) aseeauiss 
000774 210 FORMAT(LOH PRORLEM CULE 3(4€,, )1AAOe/1BH ND, OF VARTARLES 3(2me ye aMDcuNUIIA ; 
AL3e/144 NOW OF CASES 4( 2H. }15e/7249 NO, UF SELEF 1 7Ou CARUS 139/26H 9MDRINULTO H 
' A NUe OF VARTABLES ADDED 139/22M NUe OF TRNGFEN CartnS ry L3e/e2rt N  aMDEZNUIAN 

AQe OF FURMAT CARDS eH. I139//) aMNS nULAY 


c AMOS NULA? 1 
ENO AMDRIHOLR2 i 


PROGRAM LENSTH INCLUDING 1/0 BUFFERS 
008544 


{ 
| 000776 


#4DSD 
FUNCTION ASSIGNMENTS > 


STATEMENT ASSIGNMENTS is 


1 = 000136 5 = 0001] 8 - 000324 9 ~ 000332 
10 = 00062 li = (000376 le - 000404 ls = 000623 

20 = 000647 2) = 000214 2e - (00221 23 » 000564 noe 
24 = 0006490 28 = v00685 ee - y00727 31 = 060434 = 
32 = 000656 33 = QU0K37 34 *= 000765 35 - 000753 Be 
50 # 00077) 65 - 000627 yu - y00)6e 7] « 0900325 
10) = 01061 102 = 01065 lu3 - vlna? 104 = v0lnTe 
105 ~- 001100 106 = yulres lu7 = 001107 10a ~ OO1LL4 
109 - 001123 lo = vlloo?7 dll - 90)144U li? = 001143 
125 = 0012465 130 = v0034) 2v0 - v00050 20) © 0v0Ns6 

202 = 000042 203 = voonse 204 = 001043 ane = wvd0070 

‘ 206 = 000075 207 = 000NTT 209 = 40030¢ 2ln ~ 003150 

21) = 000197 215 = v00a06 242 = (00647 24a - 000593 \ 
250 = 0005¢0 283 = vourne 204 - 900706 306 = 000056 . 
309 = Qnvoboe 998 = 00224 9vy © y002?6 7777 e oodln72 


BLOCH NAMES AND LENGTHS 
- 037122 


VARIABLE ASSIGNMENTS . ie 
aj23 = 004207 8123 = ublerly cize3 = €062)1 V)?3 e 004212 ue eae 
FS = O0l2l6 FN = 004237 READ = 0 35230C0) I = 004226 ‘ not 
¢ J = 04227 JJ = 006232 K = v0423U “TA9F 4 904222 
MR = 004250 MY = 004266 NADD = 04220 NC = 004235 
NCARD = 004223 NG = 006217 NO] = u062¢ NS? = 00924}. : 
NH = 06233 NL = 004234 NINE - v0szau NNP = 00425) | 
i NP = 904216 NPV = 00623) NIAPE = 904226 Nir l = 004225 
4 | NTRAN = 008221 0 NV = 006215 NA > 003165 9 vKx (004170 
“ NY = 004236 NY) = 004256 nYy = 06265 Save =» 004214 es 
SYM = 037066C01 TOVE = £04213 x © 900HDUECD1 Ava =. 003206 
KRMAX = 004203 x4] = ovl2t6 xMIN 86) 6004244 XY = 039321001 
Y - 0042467 YMAX = yNerse YAMIN (© 9046253 é e 037}203¢Cu! 
| START OF CONSTANTS 
| 000774 
START OF TEPOMARIES 
| 001200 
START OF INDIRECTS : 
o0be2i2 : 
| UNUSED COMPILER S#ace 
' 017200 
{ 
i 
| SUBROUTINE AISTINYY NYT oAminpAmMad se SY MR QNP) RM0SD001A« 
CHrs} SUBROUTINE HAST FYB BMiTUSD IVLY 479 1964 RMOBONOLAR 
| 000011 NIMENSION XY (51017) 6A (15000) oN T4655) oXm (3) 0913) SYMP R) » amdsn0144 
KZ (15) 9 SONE (3) 9 CL AB (57) §M05000187 
' Cc amosnndlaa 
| 090611 COMMON Xo CLA oXVeSymge aMdKDNO LAY 
i 000011 TYPE INTEGER KMeiONE pUgXVod RMOaINUIaD 
600011 23 FORMATC(LM Fooly iXeAlslTAOoAl eX pF Del) AMOR DNY19] 
000011 10} FORMAT(IM SA91LG(F4.1 22X) OF 4 eo L/SIKOLO(FHL1 92K) SOK IT (atiteoen,)) BmMOnn0L92 
i 090011 1V2Z FORMATUBXe LT (OHOse, god OX elo (FG PAK) oF Gel SOKOLEIF A 1 92X)) 8M945900192 
i 000031 $000 FORMAT( BH MIN & 4FIZ,AeBUX9TH YAA & QF 12,4) AMOSNN0196 
c BATTELLE COMMENT 8405900195 
i Cc OCTAL MASKS INCREASED FOR O400 vOKU SIZE. AMOS )n0106 
000011 KM(4) # 000000777777000000U6 aMOSNNUL97 
| 000012 XM(2) ® T77777N00000000000V08 AaMog000198 
0000164 HONE (1) (61111000) ,AND.XM(2) AMOGHAULOY 
000016 BONE (2) m(+6H00N211) AND, XM()) amMDRnn0Z200 
i 000020 DOL) @(¢6QH=—"000) cANNeAM(2) amMosnnenl 
i 090022 1) (2) B(*6H000=—6) ANN eAM(1) AMDSDNUENA? 
0000246 Mal ; AMDSDCUYN: 
000025 WRITE ( He AVONIXMTNOXMAA AMDSINDZDe 
0000346 DO 50 [81935 aMOsON02NS 
0000461 SO INT(I) 890 AMOS INUSN 
000044 0 100 Kalel7 amos nen? 
000006 HO 100 Jal,50 AMDSONOZNA 
c BATIELLE COMMENT aMOS9nV209 
c HOLLEHITH BLANKS USED INSTE4U UF 0060 OCTAL 3LANKS, * BMDSDNU2IO 
000047 100 xXY¥ (Jon) @ (¢0n ) awosnodel) 
i 000057 MI NHaXx4IN/SYMS AMDKNND2 1A 
000062 TXAINSAMIN/SYMB 1 of aMOSN0219 
‘ 000066 SLAG (L) @XMIN amdennUels 
000065 00 16 18234 am04)00215 
000066 16 CLAUDE T) SCLAD (Tel) *Gy Mes aMOannveraA 
000073 WRITE ( GelV1) (CLAB(T) oABLeBeed) o (CLARY), 1eby33—2) 8M0K,900217 


0001i2 
000117 
000123 
000\e7 
000133 
000135 
000149 
000142 
000146 
090147 
000151 


000152 
000160 
000160 
000362 
000167 
000371 
000275 
000176 
000203 
600207 
000210 
000211 
090221 
000222 
000223 
000225 
000244 
000246 
000250 
000252 
000253 
000255 
000261 
000263 
900265 
000265 
000267 
000321 
000340 
00034) 


1 
I 


1 
0 


um 


DO 1 TaNnyyonyt 
HOX(T)/SSYMB=TXMEN 

INT (CK) BINT (n) o] 

TF CINT (M) =INT (kp) B49) 


MaK 

CONT TINUE 
YMOXSTNY (Mm) 
$C2d0.0 


IF (YMAK@SC) 4093493} 


sCzsCe5u.0 
60 lo 32 


WRITE ( 


6+01039SC 


FORMAT(23H TOP LEFT Hann SCaLe [> Foel) 


$C250.0/SC 
vO 6 131434 
AL@INTOL) 

Lea esce.S 
IF (L) 59605 


MBs (387T°])/641 


ITs MOD (Iee 


) 


LF (11) 92959062 


ITt¢ 


KYL MB) BLY (Lo MH) oANUG KA(AT) erReV(1T) 


Lebel 
IF (LI) e691] 
DO 19 Kelel 


KY (KoMD) BURY (Ke Md) SONU XM(ET)) CUR BUNE (TI) 


CONT TINUE, 

vO 7 K®1 50 
L23l-K« 

Ral 

RBR/SC 

Tz MOD (Ke5) 
IF (1*l) 29302 
wa lelHe) 

60 TQ 7 
aulelHe) 
WRITE ( 
WRITE ( 
RETURN 

END 


SUBPROGRAM LENGTH 


000536 
MIST 


FUNCTION ASSIGYMENTS 


STATEMENT ASSIGNMENTS 


1 

6 
123 
39 
5? 
102 


000135 2 
000246 7 
vooees 15 
voulse 31 


vo0elo 62 


000360 


BLOCK NAMES ANU LENGTAS 
037122 


Held) cow 9 (XV (19M) eMepoll) Aer 
Goll2) (Chan( Ll) odFLotheed) o(CLARIY), 1293992) 


000265 3 - vv00a63 5 » 000176 
QOU2ZET 6 - v001343 10 = 0002075 
vu0 TAKE 196 = v00040 23 = 000344 
undya? 3< = 00144 50 - 000041 
vuvell 1V0 = 000067 101 = 00035) 
000614 4000 = v0036/ 


a ae een oS eee ee NR ne ee nem am ot 


sheer 


gua n0218 
aeDdsna021¢ 
amOann02ay 
amdsNn0e2) 
3408000222 
aMDS900223 
AMDSONV22R4 
aMOSDn0Za5 
RMOSNNUCAA 
Budsonler7 
aMO50n0¢2H 
ay405900299 
8Md50N023n 
BMDRONVe 3} 
AmMDSnnle 3p 
8M05000232 
QM08900234 
8M0s000235 
am0SON0e 3K 
aMOsor02a7 
aMDRONUe IR 
aMdaNNUe 39 
BMOaHOUeEn 
BMDENNUeS) 
amMpsnn0eg? 
QMNEINOEG3 
QMDsNNUese 
amM0gn0026S 
AMDRINVC4A 
QMDENNULK? 
QMO NUZ4R 
RMNSONN 49 
AMOS nVESN 
aMDS9NU25) 
RMDENNULS2 
RMDSINUERS 
aMDSOrU25S4 
BMDEONDeS5S 
RNDRNACLSA 
amMDaNNO?S? 
AMDENNULSA 
aMDR INNS 


WVaRIABLE -SSIENMENTS 
B0NE + 090515 CLAI = 035230C01 0 ~ y0osle L - 000822 
INT = 900644 a = 000536 J - v00S26 K - 000823 
L » 900532 C] = 90082) uo - 000534 tI-4 0008525 
R = 000535 sc = 900539 51m = vs?neeCor Tkein . 900826 
] = 000520 x = 0GONN0COL Xb = 900831 AM = 000607 
xy ~ 035321001 YMax - 00527 2 = 037103C0) 
START OF CONSTANTS 
000363 
START OF TEVPOMARIES 
000426 
START OF INDIRECTS 
000440 
UNUSED COMPILEX SPacE rs 
021200 
{ 
| SUBROUTINE FLOIR (Ke ZMi No ZMAX WM TN OW MAX ONC 9 V2) aw ORDN0ZAN 
cPLuTe SUBROUTINE PLOTR FOR Havusn (MOUIFTED) JULY 205 1964 amoenn02s} 
c BATTELLE COMMEN} RAMPENNOZAD 
c MASKING VECTOR (XM) REMOVED FROM PLOTR, WOW IN FURMA, 805900263 
000012 DIMENSION VOLSN00) eV C15) oR (S101 Si oC LAR ID7) oAM( 4) oSvu(l5S) AMDSIN0264 
| Lg bF (1G) o FST IL) AMDRNNOZHS 
| 000012 INTEGER BLANKS GFoEMToTCoTPyTexyow AMOSNNNZ65 
000012 DIMENSION TR(1) eTCed) AMD&NO0eS7 
| 009012 COMMON VeCLAbexToSyveY AMDENANOZAY 
000012 COMMON /PRESET/ NCC audKINVeH? 
1 00012 DATA (NCC # 0) amdann0e7A 
| 000012 «100 FURMAT (LH OX5 (F172, 394K) oF 1203/1 IAS S(F Lec 30 5K)) ang onveT) 
000012 1NL FORMATIIH FlizetelrealslAaosnSeAl eh dee3) AMDaRNNUeT? 
| 000012 102 FOKMAT(IH 1L3A,ALe] AAG, A5eAL) auDgNNG27a 
4 000012 1000 FORMAT(IN 142101481) AMNSNNUE Ta 
car 000012 100] FOPMAT(ASXe20(5N#, 00) gIM*) QMDEnnveTs 
000012 BLANKS® (46H ) AMOSNN0eTA 
000013 IF (NCCT) 46950 04R amM05 90277 
000014 So KLS0 QMO )NUEeTH 
} 000015 GF (1) a (eoHls ) amMD5 N02 Ta 
{ 000017 GF (2) @( eoHeXr ) amDSN00eRH 
| 000020 GF (3) 8 (66HSA ) aMDanueR) 
| 000022 GF (4) 81 e6HOA ) aM aannuen? 
| 000023 GF (5) B(eOHDA ) AMDB INOeAA 
| 000025 GF (6) B( 40HOA ) QMDSONUe Aa 
H 000926 GF CFs a1 eGHIA ) BMDENNOeRS 
| 000030 GF (RB) s(+bHHX ) AMDSNANULHA 
| 00003) GF (9) 8 (GHIA ) amoannd2aR? 
| 000033 GF (10) 9 (46H 10x ) amps rndean 
: 000034 FuT( Lim (eantiza ) aAMDRINDZAD 
000036 F4T (2) RBLANKS avpsonb2do4 
| 000037 FMT(3) MHLANKS BvDEN00291 
090040 Fury (4) a(e6nS(F 126) aupsnn0zo? 
| 000041 FMT (5) a (96N39AK) 7) AMDSONOLDI 
| 000043 FUT (6) alonrlX, ’ AMDSIH02V06 
000046 FMI (BK) 6 (OHO (Flee) amognnddos 
i 000046 Ful (9) al eGhsenn) 9) RMOSNNOZ9A 
| 000067 Ful (10) aloo 12.3)) aM0490U297 
| 000051 TCa‘*lhe) RMOKNAUZIZ 
000052 TOelelee) Qmu06 700299 
000056 CALL SCALE (WMI Neoware 10 Vg JV eV MLNe YMAKO VIS) auc~ 0300 


080062 
000064 
000066 
00007Ts 
000077 
000103 
000111 
000113 
000114 
000122 
000341 
000147 
000152 
000154 
000156 
000156 
000160 
000163 
600166 
000167 
000173 
000173 
000175 
00v203 
000223 
28083) 
000237 
000245 
000247 


000306 
000314 
0003146 
000315 
000316 
000320 
000322 
000323 
000333 
00034) 
000343 
000350 
000352 
100354 
090256 
000360 
00036) 
000382 
000366 
000367 
000370 
000377 
000405 
000411 
000415 
006416 
000017 
00042) 


[ezekeke) 


223 


27 

09 

208 
i 


250 


16 
ei 


YACYMAX@YMIN 

Jevy 

TF (56510) ) 204920} 920) 
IF («L) 2200220.231 
WRITE ( 

IF (KL} 2509¢50.229 
CLAK(L)= YMIN 

N0 222 Ia¢e)) 
CLaB(I) CL ais(tel) evil : 

WRITE ( AeLOND (CL AB CT) shaletto2) o(Ches( sy ep JB201002) 
IF (KL) 2310231916 : 
TF (95) 2099221 0207 
JaJde5 

Jvrasey 

CONTINUE 

IF (KL) 2209¢26,¢e27 
FM] (3) 8GF (UY) 

Ful (7) eGF (JY) 
TIsJy 
TTRTTPViU/lve 
CLaN(1) sYMInety 

DD 223 122910 
CLab( 1) CLA (T-1) 
WRITE 

IF (KL) 227027914 
TF (SYS) 204020 992Nnk 
WRITE( él 002) 
TFi KL) 25062500226 
J2 5S @ JYT 

WRITE ( 
(TCeofeleyyT) 

TF (KL) 250,250,226 
CONTINUE 

NCC} 

1C20 

TF (NP) BUsl1e1) 

OO 1 I#) 5) 

00 ) Galet? 

MY (To J) SULANKS 
CALL SCALE (2MIND7MON 9 S000 JK OAMINOKMAKOAT J) 
XReXMAKoXMLN 

IF (WC) S2013949 

IF (NP) 80010910 
D0 9 NB) eNC 
SY4n = SYM(N) 

AD ir FREAMAKAX 

IF (ANIFFR) 105—91060,106 

XDIFFRe0.0 

YOLFF RSYMAX~Y(N) 

TFC VOLFFR) 107,108,108 

YOIFFR#U.0 

LadS2.=(50.*#XOTF FM) sane, § 
KBLOLe(LOUS*YDIFFRISYROe.9 
Kz MOU (K_6) 

Ka(K~1) /60) 

IF (mM) 21016021 
Cd.) 
Llam 

wa (me)]) #0 


4019013 


eYIy 
Hob MT) (CLAB(T) vdete}002) 9 {Clas (J) eUBL09 02) 


AelOODIITCeTalod Del CTPg (Corti gg) ) OKBL yg LVI OTD, 


BATTELLE CUMMENT 
MASK ITN PROPER CHARGCTER = SUBROUTINE ADUED 3Y aaTIrLLE TO 
REPLACE INDIANA CO@NP ROUTINE CatlLEO FORM2, 


67 


3m0e000301 
aMDannd3sor 
AMDRDNO3NA 
AaMDRINO3ne 
AMDB0NO INE 
ampnond3ne 
RMOSOOUINT 
8MD05500308 
aM059003N9 
9M960N0319 
3mM080003}) 
am0snn03)2 
QM08NN0313 
aMD5N00316 
9MD5NNN316 
aMNSNNV3I146 
3M96000317 
AMDRHHOIIA 
amMog5nn03190 
amosno03an 
QM050N032} 
3m0600032a? 
amM0son0324 
9M06900374 
aM0S0N0396 
9M9500032H 
AMOSDNG SPT 
QAMDGNN032AR 
9M05900329 
AMDENNOITA 
QM05000331 
AMDSNNN332 
8MN5190339 
3M0S9NU 336 
3495000335 
QmM05100334 
QM05000337 
AMNGOOU IIA 
8M059H0330 
aMORONOIEN 
AMD80N03461 
aM0s0n0 3407 
3M06000369 
aMOGNNE 26a 
§M08900305 
IMNRNAVIEA 
4M05590307 
am060n03eu 
aMDaonG 340 
3MDaNn034SH 
9M0690038) 
AM0R INV ISS 
aM0870U3§2 
9M0K0N0956 
9065000354 
aun 300386 
9M0GN003487 
av0a0035A 
audann03$0 
amM0annd36n 
9m06500036)1 
9mM05900362 
am0anndsar 


660923 


6354032 
000600 
000440 
600442 
060456 
06045) 
000452 
000454 
000657 
006650 
06046) 
900470 


000475 


000475 
000876 
ee0500 
000596 


000506 


000523 
000527 
000531 
00053) 
000533 
000562 
000566 
900570 
000572 
0008764 
000575 
000577 


000600 


000665 
000606 
0004607 
000611 
000613 
000617 
000621 
660623 
000646 
000650 
000082 
000656 
000676 
000704 
000705 
000706 
000707 
000716 


aAeand . 


BOO 
865 


9) 
90 


96 


97 
98 


tCcsIcel 


13 


131 
135 


CALL FORMS(RV(L eK se Lbe SY¥MH) 


CGNT TUE 

-6O TO 15 
DO be [ejol7 

KAY ( lod} BLANKS 

Lei 
50 95 Neisnc 
SY4aaS Tm (NS 
YOIFFROYMAX~Y (x) 

IF CYMRE FR) 660 2865.R05 
¥OIFFReD.0 

KNOL: LipeiaynIFERiyveeS 
we 601 (R96) 

TF (4) 90992090 
Mab 
Liem 

Kxikel)/6e] 

Ma (m=]) s 


BATTELLE COMMENT 
MASK IN PROPER CHARACTER &. “guBRQUTINE ADDED 4Y aATPeL_LE TO 
RE®Lac€ INDIANA COeNnP ROUTINE CaLLtp FORME, 


CALL FORMS IAY(LoK)s Lbs SYA) 


IF( MOD (1095) 97996097 
Walp 
GO Oo 98 
wet ; 
WRITE ¢ BelOL) Kows (KY (AON) oNa1 O17) 9M yx 
60 To 15 
MeboJK 
LLsS0em 
Tex 
gr Lee San eenerahenee 
Tad 
Diptera ranean 
wWsaTc 
K2S2 
00 31 LemMoLt 
Kako} 
T2 MOD (Les) 
IF (tl) 20302 
welp 
WRITE ( 
RLABBRLABHE XI 
walc 
GO TO 3} 
WRITE 
CONTINUE 
KLel 
60 TO 230 
NCCu0 
RETURN 

END 


GolLOL RLAB We (AT (KoN) oN) OAT) OW OL Ad 


HolLOP) WHEXT (KON) sNopol?) o@ 


SUSPROGRAM LENGTH 


a0)le6 


I eee eae 


 ampsoeGsyn 


amitann0364 
9M26000308 
2408990366 
9m03N00367 
8éM0g000348 
BMD60DNU369 


RMDG9003T) 
BMORONOIT?A 
aM060N00373 
8405000376 
amOaond376 
BudSONC3 7H 
am06900377 
&08000378 
aM089n0379 
BMDS50003An 
AMNS000381 
aedannd3ae 
9M0B900393 
am0so003R4 
SsMDe0003AR5 
amM050003R86 
au0g0003R7 
5M060¢03Aa8 
BM0SD00389 
5M05000300 
8m0600030) 
aM0S000392 
8M05000393 
8M05000396 
am0aDG03905 
eM0GN00306 
3mM05000307 
9m0506039R 
905000399 
amM06000400 


~ 8M06900601- 


8M06000402 
BM0S000403 
suMOognndens 
9050004608 
aspapngen6 
RMDRD004N7 
BMOKDNAOSOR 


- MDE DNNOeND 


§M05900410 
BMOGIN04121 
amOdann0612 
aM05000613 
AN0S000416 
9405000615 
RMORDNOSIA 
am05000617 
AMDS0004)8 
AMD8000419 
aM0S5D00420 
3M050090421 


| 
| 
| 
| 
| 


OO Ne el AAT NNN A CROCE NEA te TOA ct Rt NRE EE OT I OREO 8: AAO ORT CRETE ne A) POAC CERES TE TTT CATR RTI tn RE = pe et SO IEP CO 


PLOTR 


FUMCTIGN ASSIGNMENTS 
STATEMENT ASSIGNMENTS 
=a 


000621 
000706 
000676 
v00704 
000475 
000713 
000362 
090600 
000152 
000156 
000160 
000324 
000735 


Q01216 
00114) 
0012465 
001161 
061163 
001127 
003132 
001351] 
037303C01 
001335 


50900) TODE NEP NCoNV 9 CG 


= 000323 2 = 009656 
il = 00032¢ 13 © 060870 
16 = 000616 21 = 000617 
49 = 900350 50 = 000034 
as = 000442 90 = 900476 
97 = 000531 98 = 000533 
1902 - 000726 105 = 000361 
108 = 000370 133 = 900877 
204 = 000147 205 = 000154 
209 - 900233 2e0 = 900131 
223 - 000176 225 = 900163 
239 = 600066 231 = 00077 
965 = 00066} 1000 «= 000732 
BLOCK NAMES ANO LENGTHS 

= 037122 PRESET = 000001 
VARIABLE ASSIGNMENTS 
BLANKS = 001326 CLAB == )—- 9.352300} 
I = 0011¢2 1c = 001146 
Jy = §01236 srt = 001143 

-@ 0012606 LL = oO0llEE 
Nec = 000000C92 NP = 900000 
SYMB = P1155 T = 901131 
TT © Of1146 Vv = 900900C01 
KId ‘~ 003252 XM = 001074 
KR » 001153 xY = 035371C0} 
YJ @ 001137 ymax = 003136 
STARY OF CONSTANTS 
000712 
START OF TEMPORARIES 
ooloz2 
START OF IYNDIRECTS 
oolov2 
UNUSED COMPILER SPACE 
027400 

SUBROUTINE TRANS (NoNJONTR) 
cTRANS SUBROUTINE TRANS FUR BMD05D 
090006 OIMENSION DATA (15000) 
000006 COMMON DATA 
000006 TYPE INTEGER C1239 TODE 
000006 ASNF (KX) @ATAN (X/SORT (1, 0eA4#2)) 
¢ 

000023 Cl23e( 6H TRNGEN) 
000024 ONSNe] 
000026 MARYSO 
090030 WRITE ( 691409) 
000033 WRITE 91406: 
000037 TERRORKO 
690000 DO 1000 IsleNTR 
000046 READ{ 
00006) IF (700€"C123) 300,60300 


89 


GF 
JX 
KL 
N 

SY 
T9 
xD 
ay 
Y) 
YR 


090382 
000707 
090343 
00C4a0 
090827 
oCOT2) 
000367 
000675 
090265 
000115 
000231 
0200460 


~ 
@evuteteeaeeeaegeee 


001192 
oOLl47 
001133 
001184 
0379646C01 
003130 
001186 
001150 
003157 
001140 


“ 


TFER 
yN 
YFFR 


JULY 170 1966 


aM05000622 
9m0s000623 
9M05D90426 
9405000425 
aM0e0005a6 
aMOsto0eay 
e)s5000628 
amts000429 
BM0K000430 
aMnsna04ay 
9M06000432 
8M05000433 
8M050006346 
BMD5000435 
am0s$000436 
RM080N0437 


eee tra sateen an of enrages eee AR RR Ae en ANN a aren ent pe etn nem et maf 


000065 
000066 
000067 
000105 
000111 
cocks 
000116 
000122 
000126 
00023] 
0001346 
000136 
06014) 
000143 
000145 


009346 


000146 
000170 
000173 
000201 

- 000203 
000206 
000212 
000211 
000226 
00230 
000233 
000241 
000246 
000252 
000286 
000257 
000263 
000266 
000266 
000272 
000300 
000300 
000303 
000305 
000307 
000311 
000313 
0003146 
000325 
000325 
000327 
000340 
000340 
000302 
000345 
000387 
000357 
09036) 
000365 
000366 
200373 
000373 
090377 
090400 
on0ee2 
000404 


309 
6 


19 


24 
26 


i) 


100 
32 


NIBONS 
RETURN 
WRITE( 
MASNPNE HN 
MEGNENVENG] 
NCEMBONe} 
IF (NC® (15=NC)) 1500,1500,2 
WRITE 623406) 
80 YO 1000 
IF(NC#11) + 3, 3 
Kalo 
MDENSKON 
00 200 JemBeMC 
MAamAe) 
MOSMD+) 
CONTINUE = 
BO TO (10220030 060.50060070,800909100110212001300140) ONC 
IF (DATA(J) 9903208 
DATA(MA)SSORT cOaTacu); 
60 To 200 
IF (DATA(J) 99611912 
DATA(MA) 81.0 
60 TO g00 
DATA(MA)SSORT (OATA(J)) eSORT (Oat A(U) 0100) 
60 TO 200 
JE (OATA(J)) 99990016 
DATA (HA) BALOG (DATA (YJ) ) #0 046302944619 
60 TO 200 
DATA(MA) SEXP (9aTA(Y)) 
60 TO 2c0 
TF (OATA(J)) 99932017 
IF (DATA (J) ©] 00) 18919099 
DATA(MA) 03.1618926536/2,0 
60 YO 200 
ABSQRT (DATA(V)? 
DATA (MA) BASAF (a) 
60 TO 200 
A@DATA(U)/0N 
eho} da 2 
F(A)99» 23,24 
ae 
DATA(MA) 80,0 
@0 To 200 
DATA(MA) BASNF (SQRT (8)) 
60 TO 200 
1F (8) 99,28,29 
DATA (MA) @ASNF (SORT (A)) 
60 TO 200 
AsSQRT (A) 
QsSQRT (8) 
DATA (MA) BASNF (A) CASNF(B) 
GO TO 200 
IF (DATA (J) 231099931 
DATA(MA) @3,0/DaTA(y) 
60 TO 200 
DATA(MA) BDATA( J) oCO 
GO To 200 
DATA (HA) MDATA( 3) 8CO 
eo To 200 
IF (DATACI)) 33932033 
DATA(MA) 80,0 
ec To 209 


601402) 1 9NE ONC EN CO 


QNDS§000498 


en0s000439 
9m0g0004660 
omdg0004e1 
aMognn04e2 
9m05000443 
am0s009444 
am050004465 
8M050004466 
9m08000467 
emd3000468 
84050004490 
8n050904680 
6m0s00048) 
andnng04g2 
am0s000683 
Bmdgnndsse 
am0s00045s 
aM0en004K6 
em0s000487 
oM0s0004s¢8 
aM0g000489 
9M08000460 
M08 00046; 
908000462 
9M05000463 
8M05000464 
9M05090465 
8M05000466 
9M05000467 
pmdenn046a 
0m05000669 
6M05000470 
9m08000471 
9M05000672 
8M08009473 
OmM050004674 
808000678 
Bmb8000476 
9m08000677 
0405000478 
608000479 
em05000480 
8M0800048)} 
8u0so004a2 
9m080004683 
8m0S000084 
0M05000488 
aM0S000486 
BM0S000487 
oM0s000688 
640500049 
6m08000490 
0408000691 
9m050004902 
8405000493 
9408000494 
8405000498 
BM05000496 
6405000497 
9M05000406 
9405000499 


to et see erent etnies! 


a ggaremcnctancnayel 


wee 


000495 
000413 
000416 
000422 
000422 
000430 
960436 
600435 
000436 
600440 
000445 
000446 
020450 
00065) 
000452 
000460 
000466 
000673 
000474 
000477 
000500 
000504 


000505 
000505 
000505 


000505 
000505 
000505 


000505 
000505 


0005805 


33 DATA(MA) BDATA( J} ©*C6 
eo To 200 

116 
@0 To 


@0 16 


90 To 200 


140 IF CCATACMO) 11570990157 


157 DATA(MA) BDATAY 5) (DATA (MD) 


60 TO 200 
99 TF (MARY) 63066044 
46 MARY2~999 


BATAMar SDATAC ICDS COND) 
DAYS A SOATA(JT=ORTS (HO) 
GATACMA) @DATA( J) @DATACHDD 


TERROR®=999 
WRITE ( 6e1400)1 
43 WRITE( $91405)0 
200 CONTINUE 
MARYa0 


1000 CONTINUE 
IF (TERROR) 4251, 3103112 


42 WRITE ( 601401) 
1321 RETURN 
900 FORMAT (AGeI3Seloel39F 600) 
Cc 
1600 FORMAT (4SHOCARD NEW 
YVARTIABLE Code vaRta) 


1601 FORMAT(7@H VALUES OF VARIASLES Of WHICH AN ERROR WAe FOUND DURING 
LJTRANS@=GENERATION WILL /77H STILL BE INGLUDED IN THE BRAPHS, HUWEVE 
2Re THESE GRAPHS MAY BE MEANINGLESS 


TRANS 
OR CONSTanT) 


3€ TRANSFORMED and OTHERS NOT.) 


1602 FORMAT (2H 


12918¢21906X oF 10,5) 


ORIGe 


/S4H STNCTE SOMe VALUES WILL 8 


1603 FORMAT (IHO6Xe23HTRANS GENERATOR CARD (S)) 


1606 FORMAT(SSHOTHE INSTRUCTIONS INDICATED ON TRaNS GENE@aTOR CaRD NO,1 
1201X%)3HRE=/60H SULTED IN THE VIZLALION OF & RESTRICTION FOR THIS T 
GRA RSF OR@/aoH MATION, THe VIOLATION OCCURRED FOR THE 7TEMS LISTED M 


1605 FORMAT(10OM ITEM NO, 15) 


14606 FORMAT(LOTHOTRANSGENERATION CODE ON CARD LISTED ABOve IS INCORRECT 


Wine Van(6)/45H NO, 


Ke PROGRAM WILL PROCEED WITHOUT THIS TraNSGENERATION,) 


€ND 


SUBPROGRAM LENOTH 


000732 
TRANS 


FUNCTION ASSIGNMENTS 


ASNF 


STATEMENT 


000010 

assraent 

000131 - 006134 
000067 : = 000173 
000211 14 - 900233 
000263 20 = 000203 
000323) 2? = 000314 
000239 31 = 90036) 
000264 42 = 000800 
000254 60 = 000300 
000373 99 = 000646 
0004822 190 = 0004630 
0004666 300 #» 900065 
000504 1900 ==) (000824 
000574 1904 © 000601 
000122 


aeaeeceaeuace ae @ aan 
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000141 
000170 
000257 
000307 
000327 
000402 
000460 
000357 
000400 
000436 
00082) 
000840 
000626 


000146 
000206 
000266 
000325 
009340 
009405 
000450 
600366 
000614 
00044" 
0004674 
000870 
060631 


om05500500 
3”08000501 
0405000502 
6408000503 
anM05000504 
8405000505 
oM0S0005S06 
6M05000507 
em0so00808 
9m05000509 
embso000Si0 
6m05000511 
6m08000812 
am05000513 
81050008146 
am9s000515 
@m05000516 
8405000317 
8ud35000S1a 
8408500519 
9m05000820 
8M0800052) 
8408000522 
8™05000523 
@m05000526 
amD8000525 
9u05000526 
3n05000827 
BMDgN00S2a 
aM08p00829 
9™05900530 
2405000531 
8405000532 
9M05000533 
8M0S0005346 
9405009535 
8M05000536 
9408000597 
905000538 
RM05000530 
8405000540 
840800054} 


@LOCK WANES ANO LENGTHS in 8 a 


FUNCTION ASSIGNMENTS 


STATEMENT ASSIGNMENTS 
10 = 000004 20 = 000013 50 * od00le 4000 = 000022 


BLOCK NAMES AND LENGTHS 
VARIABLE ASSIGNMENTS 


START Of CONSTANTS 
000020 


START OF TEMPORARIES 
000034 


© 03523¢ 
VARTARCE ASSIBNMENTS 
= o6CT3¢ 8 = 000731 co © 060721 C1232 w 000710 | 
oatA = 000900C01 1 = 000715 TERROR = 000714 J © 006727 . 
K = 000723 wa = 000722 MARY = 00713 "6 =~ 000723 Rae 
mC = 000726 “od = 000726 NC = Qn07l7 NE e 000716 ven 
NY = 000726 ON = 900712 TOpDE == Ss O0TIA pots 
START OF CONSTANTS eae as 
200507 ae 
START OF TEMPORARIES Cate 
000646 aoe 
START OF INDIRECTS eae 
| 000792 ay 
UNUSED COMPILER SPace Poo 
020400 eal 
| SUBROUTINE VFCHCK (NVE) E e40S000S42 aren 
CVPCHCK SUBROUTINE TO CHECK FOR PROPER NUMBER OF VaRTARLE FORMAT CRos 998000849 
| 090003 id Soyer Lei Oize SMDSHNOS 44 oma 
000006 19) WRITES 624000) 8MD80N0S4g bo 
900080 nVFel 8M05000566 : 
000012 50 = RETURN 8M08000547 Bee 
| : supsp905es a 
F(NVF210)50950010 M0800 . 
| 000013 Re 1 710)50250% euoebaeess ie 
000616 #000 FORWAT (140237 LHNUHAER OF VARIABLE FORMAT Caa0S IncARRECTLY SPECIE 9M08D0055} te 
| KJEb» ASSUMED To BE 1.) RM0SD00582 pas 
000016 END 9M08000589 aera t 
| SUBPROGRAM LENGTH ail 
; 000036 eed 
| VPCHCK a “ 
| : 
| ; 
i aoe oe 
a 
' . 
} 
| 
| 


START OF INDIRECTS 
000036 


UNUSED COMPILER SPace 
022600 
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1 TT TREE sean tn ee aOR Yet ARNORNERRE th fenabyninen tea wee ot cere 


SUBROUTINE SCALE CYMING YMAAGVINT IGT eo TYMINGTYWAK. VY Jy 808000556 
CECALE SUBROUTINE SCALE FoR Sud PLOTR AUGUST 18s 1956 8MDSD00585 i 
000012 DIMENSION C(1@) BMOSO098SS 
000012 Ccho# 3.9 9M0500 0587 
000013 C¢(2191,5 au05D00588 
900014 C(3)82,0 aM0500088¢ 
000016 Ci4) 23.0 8408000860 
009017 C(5) 86.0 3405050546) 
090021 €(6)e5e0 asdSO00S62 
000922 C{7) 8765 808000563 
060026 C14) 81000 BM9S000346 
| 000025 FESTa, $#(200(@29)) 8408000565 
000033 5p YAsYMAXeYMIN AMD 000566 
: 000034 TTeYR/YINT 8é05000367 
000036 JBALOG(TT) 064302934819 3™05000868 
000062 Ea) 9,0¢e) 8408000569 
009046 TTstTt/€ 8M05000570 
000047 180 805000571 
600050 IF (TT#220)205520).201 RMDS000S72 
i 000056 205 TTett#10.6 eM0S006S73 
i 00060 £96100 BM08 0005874 
000062 201 Istel 8m05000875 
000066 IF (Gal) leee2 8M0S000S76 
000066 1 €5E®1000 3M08900877 | 
| 606070 Isl RM0500057A 
000073 2 LFCTT$C(1))2330200-20) 8M0%000579 
000075 233 YIUsC(1) *€ am0so00Sa0 H 
00010) 60 TO 203 9405000581 
000102 202 YSYMIN/C(Y) aM050005A7 
000205 Jay 9M05000563 
| 000107 Tay AM05000586 | 
000110 IF(0,0001-ABS (Tov) ) 204923390233 RM05000568 
000115 206 YVidec( let) E aMds000584 
; 000120 203 Xa((YMAKCYMIN)/YI JeVINT 172006,00001 \ amM)5000887 
000225 ax RM05000S88 
000127 TECK) 23592400200 am0a0005a¢e ! 
00013) 235 Yak RMDS000590 
000132 If (XY) 23692490236 RM0500059} 
000135 236 Kek=l 8M05000592 
000137 260 TYMINSK am0s000594 | 
000140 TYMTNGYLU®TYMIN aM0§000896 | 
000142 TYMAXaTYMINGVINT@YI YU 8M05000508 
000184 IF (YMAX@TYHAX@TEST) 100100¢01 AMDSO C0594 
0001467 lo TTSyYINT/106 805000597 
000150 JYsTT*.000001 RM08000S89R 
000153 YIJeYINT@(Y¥ 1371000) am05000899 \ 
000155 JBTYMINS YIU aM0K000600 
000157 IF (K)2420ee1,261 aM06000601 
60016) 242 JaJ@l am05900602 
000163 24) Jade eJYaK 8M08000603 
"90167 JY¥ey. aMd50N0606 | 
00c\70 RETURN 9M05N00605 
000176 END aM0SD006N6 | 
SUBPROGRAM LENGTH 
000206 
SCALE 


FUNCTIUN ASSIGNMENTS 


STATEMENT ASSIGNMENTS 
1 = 000066 2 


00007) 10 000147 50 « 000033 


201 ~ 009062 202 = 000102 203 - 000120 206 - 000315 
205 = 000056 233 # 000075 235 = 600131 236 = 000135 
260 = 000137 2¢1 = 000163 e*2 - 000141 


BLOCK NAMES AND LENGTHS 
VAKIABLE ASSIGNMENTS 
¢ 


= 000720 e = 000236 I = 009237 J - 000238 

K = 000243 T = 90024) Test = 90023e tT - 000234 

x = 000202 Y = 00260 viy = vo0Go00 Ya » 000333 

START OF CONSTANTS | 

000172 
| STaRT OF TEMPORARIES 

000211 
| STaRT OF INOLRECTS 
i 000217 | 
UNUSED COMPILER SPacEe 
021700 ag 
i 
| 
SUBROUTINE FILL(SY4s NC) A¥05000607 | 
c AMDgD006NA 
C BATTELLE SUSROUTINE AMD50006N9 | 

c RMDRNDNbIA i 

Cc THIS ROUTINE FILLS THE COMPUTER SURD WITH THE INPUT GYMBOL SYM(N)  9mM08000611 
} Cc ADDED TO INDIANA CO=0F VERSION IN ORDER TO COM>; [MET ADDED aM05000612 
c BATTELLE SUBROUTINE CALLED FORM36 AMD5000612 
l Pe 7. AMDSNNI616 
¢ THIS VERSION IS FOR & XeRYTE WORUs WHERE X 10, BMD&000615 | — 
c aM05900614 — 
600005 DIMENSION SYM(15) gM0R000817 .g 
¢ QM0500061A 
090005 00 10 NwlyNC BuDRDN0619 
000006 10 ENCODE (1091U009SYM(N)) (SYM(N)» £81039) AM0g000620 ; 
000031 RETURN am0$00062) i 
c AM08000622 
00003) 10006 FORM.T(JOR1) AMOS 000629 
000031 END AMDSD00626 es 
SUBPROGRAM LENGTH 
000043 
i FILL 


FUNCTION ASSIGNMENTS 


STATEMENT ASSIGNMENTS 
1000 = 000034 


BLOCK NAMES AND LENGTHS 


VARIABLE ASSIGNMENTS 
1 = 000042 N = 000061 


| 
START OF CONSTANTS . 
000033 —— 
START OF TEMPORARIES 
74 
’ 


000036 


i START OF INDIRECTS 
00004) 


UNUSED COMPILER SPace 
022500 


ccc Heading aveste 


000006 
060006 
000006 


099006 


000006 


000006 


000010 
000012 


000018 
000020 


060022 
000027 


000030 
000032 
000035 


000060 
000005 


000046 


AAINNAANANAN AANA ANAAAAMAAMAANnNA 


aa aan 


an 


\o 


24 


SUBROUTINE FURM3(GRaPH, LOL, SYMaUL) am08000625 
AmMDgda0e24 

BATTELLE SUdROUTING 5M08)00627 
8MDRDNG62A 

THIS SUBROUTINE PERFORMS ALL THE MASKING OPFRATTUNS FOR PLOTRe AMDRNNNE29 
amdsno0e3n 

GRAPH © THE COMPUTER ®ORD AT THE CUORDINATES 19 SF , NADEU. 3M05000631 
LOC @ THE BYTE LOCATION WLIHIN IHE COMPUTE? 408n Foo LOADING. 3M0R0N0632 
SY480L = THe INPUF CHARACTER TO 32 LIaneuU. 481 Flou Uwe 9405000639 
SIA LEFT 8CU CHARACTEKS OF THe CUMPUTER WORD. AMDS000634 
QM05000635 

THIS ROUTINE AND THE ENTIRE PHOGHAN TREAIS T4E PLOTTING WORD AMIS DN0624 
AS IF IT WERE ONLY & ALPHANUMERTC CHARACTERS LONG. | eR T-JUSTIP AEN, aM06000637 
3mM060n063A 

THIS ROUTINE LOADS THE INPUT SYuauL aT THE | ICATION LOC WITHIN §m05000639 
THE COMPUTER WOKD GRAPH, AMDKONNGEN 
IF AN JIEM ALREADY FXISTS AT A LUCATION TU Ef LnanEn, 8MD50N0646} 
THE ROUTINE WILL LOAD aX NUMBER UN LETTER INDICATING THE NUMBER aMDGD00642 
OF {TEMS AT THal LOralION. MEFER 10 TanlLe FIX MULTTOLE@L lem 905000663 
ALPHANUMERICSe Tre NUMRER OF [TemSe Ny 1S F2Uar TO TRE aMORDN0644 
CHARACTER STORED AT TABLE(N@]). aM0G9N0645 
NOT ALL CALLING PROGRAMS USE THIS SULTTFLE@T TE FEATURE, BUT 2M057006446 
IT 4AS HEEN INCLUVEN FOR THE Sane UF STANDAQ)EZA1{Ou, 3698900647 
BMONDN0G4H 

DIMENSION TABLE (15), MASK(0) RMIEDNVOAG 
INTEGER TABLE® GRAPH» SYMBUL» dL ANAS QaMdaNn0eSe 
DATA (BLANKS#0H ) 8MD50N068) 
aMDaDAD6R2 

DATA L TABLE s6H222222 » 6933353 1 OHS 44446 GHS5 55559 6HOOKKOOOOKTTTI SS, B8MDG0n06S3 
\ SHBBABBS » 6H999999s DHARARAR, GH ISB KACCCCCs OHDDDUUU, AMDRDF06S4 
2 OHEEEEE Rs GNFFFFPE GHG 2G00S, BHNNHHIM yg CATT TTI yOMs/////) RMI60NbO55 
RM0S)N0456 

DATA (MASK & T7YUDNDNNOVOQVOVOUOONVEs OOTTIOONIVONVOONHOVOSs aMD05 000687 
1 QNOO?THULONOOVOOUVOULHs DOQH0VT7NDNNVOONNOVOHs RMDS50N065A 
rd VOUOONOVTTOODOUVVVOUOUHs DONNVODDDO7T700NANAUUA) am0c0n0659 
amDe)0066n 

8M0500066} 

MASK OUT LOCATION UNDER INVEST Liat lon, AMD6500662 
ITEM ® GHAPHe ANID oMASK (LOC) 2MD5000663 
AMORONDOGe 

TES! TO SEE IF IMERF ARF ALREADY 19 OR MURE LIENS ay LOCATION 9M0K59N0666 
LOGIC 8 TABLE (18) oe AnD eMASK (LOC) 9M0S5 000664 
IF CITEMLEQ. LOGIC) RETURN amM05nUOG7 
aMDS50N066H 

TEST TO SEE IF LAM (Tee NO PREVIUUS TTEMS al 1 OCATTON) « RMDRNANGAD 
LOGIC @ BLANKS, ANU MASK (LOC) amOdgnndorn 
IF (CITEM.NEeLUGIC) 69 IN 10 am05)0067) 
amMdan0e72 

IF WLANK AT LOCATIONs MASK IN SY4DUL. 9495900672 
GRAPH B (GRAPH LANDS (NUT eMASK (LUC))) ORs (SYMSOL LAND, MASK (LOC) ) aMO50N0674 
WE TURN AMOSNHVOTS 
aMDNSNNOOTA 

SCAN THRU MULTTPLE@TIEM TAOLE Fox MAlcH. gu0sdn0o77 
O00 20 Talel! 3MDGqNNGOTH 
LOGIC & TAHLE (17) ANNs MASK (LOC) QM0SOn067R 
TF CI TEMeNE LOGIC) 69 TO 20 aMOsnndoan 
aMDann0oarr 

IF MULT IPLE*SYMMOL TS FOUNUs MASS LN NEX! Sy4d9) e aMDEINDOKD 
GRAPH & (GRAPH ANU, (NOT MASK (LIC) )) ORS (TANLE (1911 BND MASKI(LUC)) aM08NNUHRG 
RE TUQN AMDRINDOHS 
amOanndors 


CONTINUE 


aMDaNN GORA 


i 


z : the a Renee Pa : SUPE ey rhage ae 
; \ septate Cecane See sits GE ah SS AE NS oe ee eons ak Tey cette et pce na satin uh Shee icen Sia tage ie, ai cee ates eS Soe eat edie aha seat Ce nee ie snap IE ate oe he ee 


FUNCTION ASSIGNMENTS 


Cc IF NO MATCH IS FOUNDs THEN MASK Iw a 2, amDs0n0687 

000050 GRAPH ® (GRAPH. AND, (oNGTeMASKILC))) 2ORs (TAILE (1) 2 AWM MASK (LOG) ) 9mM0500.06A8R 
000056 RETURN amdg0NOHRG 
000056 end awdsnn0e9n 
SUBPROGRAM LENGTH 
060131 

i 

} 

| 

{ 

| 

| tUKM3 

{ 

| 

| 


STATEMENT ASSIGNMENTS 
19 = 000030 20 = 000Na6 


BLOCK NAMES AND LENGTHS 

VARIABLE ASSIGNMENTS 

BLANKS = 000125 I = 000330 ITEM = 006326 bOnlC © 000127 
MASK © )§=—(000117 TAGDLE = 000075 


START OF CONSTANTS 
000060 


START OF TEWPORARIES 
00006) 


START OF INDIRECTS 
09007) 


UNUSED COMPILE S$Pace 
022300 
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600003 


000003 
006003 


000003 
"20003 
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C4HwW0eD CURMFLATTUN wii SRANSOENF 4a! fuw "MOVEMBER JSe 19466 
¢ O400 CUNyvERSTON arliv +ENenN 
c HE OE TR SCIENCES SUMPUTING BaciLrtiy, icha 


100 


oe 13) 


a9 


TYPE INTEGER WK ASK HI2 Fede Z OCI CSOUISI9O] TIVE CODE REL eOSeF MT 901 
is 

TYPE INTEGER PO79PO] oF 02 PUR 

DIMENSION AMAX (150) AM] (1590) 9YV(E9) 9S¥YM( 120) ONE (I An) oSK6150) 
4 9SA2(190) oSAY (1356135) Al a(1So?) oFmt (yeu) 

NEMENSTON NOW(150), JUMPC1SU) gNA(L5U) 9AN(I50) 

DIMENSION CON (36) eNSub (36) OHEL (30) 909 (96) gp INIEK (36) gwk (A) 
OIMENSLON fhaTa cle) 

COMMON SXY 

COMMON TDATA 

UVIMENSTUN C1255) 90150) 


FORMAT ISGHIeMUNED COKKEL ATLOW WI'M THANSUENF ATION © VERSION UF 
AISHNOVe 139 L9he / 

1e04 HEALTH SCIFWCES CUMPUTING FacdbiTy,ULLas/ 

C144 PXOGLEM CONE A&e/ 

3214 NUMBER UF vaRlARLES 1397 

©1744 NUMBER UF CASES 59/7) 


PC72(+24NU) 
wkK (]) aleonGl 
wk (2) 2 (e6nbe 


wh (3) alegre! Rios 
WK (4) 46HLE ye 
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WK (6) mC + 6rNE R 


wh (7) a lesnon 

wK (6) S(e+6HAY 
ASK8 (6 nt ) 
Ci e@se (SH TRNUEN) 
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H12Jse (¢6HPRONLM) 
Al232(¢modh IN (SH) 
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NTaPEZS 

Tle 

1Tie} (O2TGTr aL) 
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VE CODE (Sue TUee TUATA) TUM e CU eNVANDNSAMSNSEL sNANUDNANNL PW] ePwaePQ 
L9eNiGeMlaAPeEskvA 

HE WIND SE 

LF UINNESAL SS) (008M) 100 

IF (TOVE#31 25) 901097 3ISIDULG 

LF (MTAPELEN et) 7414/42 

MTAPESD & 90 10 773 

[F (MTAPE NES) GO TD 773 

PMHINT 5 

FORMAT C/7/////7% YOU HAVE aSS[G VEU LOGILAL WUweEa « 10 YOyR vATa 
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000155 
000157 
000157 
000164 
000167 
000176 
000201 
000203 
000215 
000220 
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00u226 
090232 
000234 
000251 
000253 
00027) 
000274 
000301 
000305 
000307 
000312 
000313 
000316 
000321 
000342 
000346 
000350 
00035) 
000365 
000370 
000373 
000375 
000490] 
000402 
0ou4e)) 
000411 
000413 
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000423 
000425 
000426 
000436 
000436 
00044) 
000445 
OCuse? 
000459 
00045) 
000452 
000453 
00046546 
000455 
000456 
000457 
000460 
00046} 
090402 
000463 
000464 
000465 
000466 
000467 
00047) 


stop 

CONTINUE 

LF UINVAR@E PO (1L36eNVOR)) 9O4UsINGUFLO0OR 

EF (nSAme]) 9001, 9007-1007 

TE C(NVAR+NAUD@17) 8 0} FOONVAHENAUD)) FU0299162, 705 

{F (JARS iNRUUL) @9) 766979699003 

CALL VECHCK (KVR) 

WRITE ( He lV) COUE eV 4 teNS aM 

IF (NTS*ESU) 1H03e 110394004 

IF (NTG) 90059401 6402 

WRITE ¢ 64403) 

WRITE ¢ 48) 

00 707 TeleN{G 

RE AU { 59406) TIDESNEW( L) eYUMP(T) oNACL) SONG TY 
IF (TANE=*C123) 9ndbesnSegnne 

WRT IE ¢ HeGO 7 To NEW( TY SIIMP(T) ONAC I) ew {T) 
IF (UP 009991) 2000070792009 

LF CJUIMP CE) @ OL Te guMP tT) ) 2005020039 fU7 

WRITE ( &,4U01) 

JUMP (J) 899 

CONTINGE 

TF (NHOVL) &119461296)1 

KR STAHS (nNKHOOL) &4 

WATTLE ( &412) 

HE aU ( 50414) (SUS CL) e4el (1) eCON(T) 00571) ore) oKK) 
wRIIE ( 64415) 

0 16 LaleKrK 

KKisl 

WRITE ( HeGL7INSUH CL) ott (1) eCON(T) 99E CT) 


LF (ASK *0P(7)) 416919969416 
CONT [NUE 

09 438 (aleXnl 

INDEX (1) 80 

10 437 Lwl earl 
DECUDE (Le 704% sCON([)) PUS 
FORMAT (AL) 

IF (AL@PUS) Tlde TLL e710 

VE CONE (59712 eCON( 1) ) CONC) 
FORMAT (2X eF 300) 

{NOEX (Lal 

60 TO #37 
VECODE (oo 7) 39CON(E 2) COW(T) 
FoAmAT(FO.0) 

CON INUJE 

IF INIGYNBOOL) 423946169425 
LF (NTG) 400794394826 

IF (NBOVL ) 42294219422 
JESUSa1 

60 To T 

JESUSH? 

NOS) 

ao 10 7 

JESUSa3 

bo to 7 

JESUS ae 

NUOJan 

G) To 7 

JESUSs> 

NO380 

mzQ 
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00053) 
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000574 
000604 
200606 
000606 
000622 
00u626 
000625 
000626 
000642 
000044 


00066) 
00067} 


000706 
000716 
oov720 
000733 
C00735 
000737 
00074) 
000742 
000743 
000746 
000750 
000752 
O0UTSs 
000T6) 
000762 
000765 
000771 
000775 
000777 
001003 
001005 
001006 
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° 
o0l026 
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001026 
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001067 
001054 
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7 
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43) 
&<B8 
09 
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433 
436 
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101 
43 
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4005 


Lous 
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onee 


eane 
eoen 


i901 
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sus 


AMIN(T) 810,019 
AMAA({) woamIN(T) 
S4() 00.0 

SX2(1) 8006 

DO & J8LleNVarl 
SXY(C IgV) 2020 
KLSvU 

fide 

KY REK VRS 

READ ( 

00 600 I1=}e9NWSam 
RE AL ( MPAPE Fm) (ATACL) gp LE) gNVARS 

GO TO (00764274428 5625,650) ESS 

CALL COUL (NSUG pREL eCUN DOP o ANDERS VATASNTEST 9 1K 2KK19A1 93 9BL299DLEIAN 
11 APE) 

GO 10 (6009431 6499e701) oNTEST 

NOSanvdse] 

GO TO %u7 ; 

CALL TRNGEN (DATA sNVAWONTUSINSAMOLER | 9LCASE ONE AD SIMD ona gBNOMe IL? 

IF (LCKS5E) 40499607 9607 

LCASEs0 

GO TO 6uU0 

CALL TRYGEN (VATA ONVOH ONT ONSOMOLEP TyLCASL ONE Ae SIMD eNAgSNoMe LP) 
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VAHIABLE ASSIGNMENTS 
AMAX = 002612 adn = glebeyu ASK = 4b26TS al = Gvelo2 
Al23> =) «(UDO f6 aN = y0527) H223 - d020T7F C = 005673 
CODE + on210a CON = y05517 Cles - vu2iov vara 2 00e,el | 
pie3 = us2yua Fai = vl221s FN = 00655¢ A © 606845 
fala] = 006559 1 = ynokie TUatA = veReGlCO 17 © 006846 
INDEX = 005627 Th = y068ei J = 1065463 JEcUS » 606835 
Jump > ~p0ebls Ke = yl6533 KAl = 06536 RL = 006544 
KVR © wu96S3} LOASE = uv6SeU Leet - vunses ~ 006837 
MTAPE © 496530 na = ytdna3 NADH = 06525 NINJL © 006826 
wEw = 004367 NI = vl6ass6 Nw = 906555 NOQu = 006536 
NPAGE = 16553 NS =~ yu3275 N>a4 - 006524 NSeu - 906824 
N$SUB = ns503 NiaPE = yl6520 atest = yoeser NTr. ~- 0068527 
NYAR * 0905¢2 NVAR] = uUloRee? or - 002151 Por = 9026406 ; 
PQ] » 002607 Pue = vyi2elu Pus © udPele PUc - 002405 
9 = vudbele Qe = Yh6sst REL. - vu2)ns SX ~ 0034665 
SKY - OV0uv0ubECO! Ske =~ vt3713 SY¥u = yv0310>5 fone = 00€103 
aK = UN2065 v¥ = Hwusnes 
START OF CONSTANTS 
001611 
STAKT OF TE WPOxARIES 
002011 
START OF INDIRECTS 
00205) 
UNUSED COMPILE SPace . 
014700 ; 
SUGHOUTINE COOL (NSUR DHE, sCUNGUP s INUVEXSNATASNTESTOWK KK LGALAS~012%_ «IMNPINUIAS 
iDlLedesNf are) AMDPINVANG 
ccoe, SUSROCUTINEG CUOl FUK seid QMO AN NUan7? 
c ReWKITIEN GY fy BUTS J=25664 AMOZINU SOR 
Cc aMdANHV3I0G 
Cc DEFINITION OF NITES: aMor9nUa1O | 
C NTEST#)] IF CASE FAILS KNOLGAN Ted am)2200311 | 
Cc NYESTw2 IF CASF SATTSFYES SNULEAY TEST ampzn0n0312 
c NTESTa3 IF NEw PROS, EM CARU IS HetecliFn am02900313 | 
c NTESTs@ TF FINTSH CAaU TS VETECTOU amMOr00314 i , 
¢ AM)? 500316 | : 
000020 TYPE INTEGENM WK eOP REL es lode bl2reui2aaX RMIANNUITA | 
000020 UTYENSION NSUb (36) EL. (95) 9 CONE 36) 9079 (36) pT VIEK (3H) yak ( &) DAIA(TS =«8MD2ONDII7 | 
Lode SXYV(13501395) eIN¢ 37), EDATAC LS) AMI2N003) 4H 
Cc AMUPDNOITE 
000020 COMMON SXV_1VATA amMD200NI2N 
000020 vO 100 im} KK} 9éM02N0032) 
00002) 1SenSud(T) 9M)290U322 
Cc Qu 0700323 
c EXAMINE BUQLEAN EXPRESSION FOR Glee GEo Lis LEe ely we . AMO INU 3246 
c 9mM029003A46 
000023 OO 55 JB] 96 AMDPONUIAH 
000025 1F (REL (CE) eek (0) ) 55520955 BMD 20N0377 
00003) 5S CONTINUE AMDANAV IA! 
000033 GO To 3)] aMDA NU I29 
000034 26 IF CINDEACI)) 27027928 aMD>9NU SIN ’ 
900037 27 CCECON(L) a40290033) ay 
090042 GO 10 29 AMD2DNU SIZ 
0000462 28 naecon(]) QmMfunv3 33 
000045 CCBUATA(K) . RMD717N03 Ya | 
000067 29 wevala(ts)eCC AaAmD2900335 | 
. : 
1¢ 


000082 
60006) 
000071 
290103 
e00115 
000117 
000121 
000123 
000125 
060127 
000131 
000133 
000135 
090137 
00014) 
00043 
000146 
060147 
000147 
000151 
0001546 
000155 
000157 


000160 
000162 
000166 
000170 
000173 
000173 
000177 
000203 
000206 
000210 
000212 
000212 
000214 


000214 
0002)7 
00ual17 
000222 
000230 
000237 
00024) 
000253 
000261 
000263 
000266 
000270 
000270 
000272 
ooveT2 
000276 
000275 


gAaananan 


Sol 


19] 
199 
222 
223 
200 
900 
32) 


320 


341 


301 
313 


302 
31? 


3u5 
307 
306 


309 
306 


LF CEES eE Dc He) eNews (ME TA(TS) oF Mell I 229G,27 300 


amaGs (6/CC} 248 (-35) 
GO TOL 293049Seb) oJ 
60 TW LLel2sl3or4etSsiny os 
IF (#)S50eS5U,11 

LF (AR) S0950 0c) 

LF (4) 20920912 

IF (3) 50020920 

TF (8) 50084913 

IF (6) 20-50_50 

IF (a) 200e20s1h4 

IF (t:) 20920050 

IF (A) 20920,15 

IF (6)50e20 050) 

IF (a) 50950515 

IF (i) 20s50—9¢0 
IN(T) =] 

GOD TO £00 

TN(T} 8d 

CONT TINUE 

NTESTsIN(1) 

KKakKiel 

LF (KK) S00250usSu) 


ExXan{(NE AQOLEAN OPFRAIDR FUR ORS AN 


HO 200 L#¥leak 

IF COP (CL) wk (7) ) 22291940222 
IF (WTEST) L¥Y¥9 t 99937) 
NTESTRBIN( +L) 

GO TO 200 

LF (OP CT) mw (3) ) 3016723930) 
NTESTRNTEST#IN(T +L) 

CONT TINUE 

JF (NTEST) 320 2320932) 
NTESIz2 

GO To 333 

NTESTs] 

GO 10 333 


ERROR LUUK FUR NEAT PROKLEM OR F LINES CARD 


AL2I26MF INTSH 
H12326HPRUBLM 
VL23ISS4PLUTSL 


K=MEL (1) 

6 TO 313 

x2QP(1) 

WRITE ( SeclOOIA 
LF (NTAPF HS) 30995026304 
JaNnTAPE 

RE AU(J91000) (ILNalTA(h) eKaeleda) 
IF (LD 4TAeU1 23) WSe 3b 29309 
LF (POAT a=61 3) 0793669 WE? 
IF (POAT AeA eS) Bee 390312 
NTEST83 

60 TO 333 

NTEST&6 

60 10 333 

REAEND NIAPE 
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60 To 312 


RMDADNVI 36 
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RMDP ONUI39q 
RMDADOUI466 
3M020036) 
am02009342 
ampennd ss: 
amOrnn0 see 
AMDr0GU 345 
RMD?DNG 344 
RMOPNNU 347 
BMD200U3aa 
aMD2000369 
AMDADOUSSA 
AMO?2ONV351 
AMDPD UIE? 
aM0?70NU353. 
aMO20nUS56 
852000355 
AMDZNNUISH 
aMD2NNU3S7 
AMDANNVISK 
8M02100359 
amMD2nnb364 
2402000361 
3¢07D00367 
RMOZADNVIB® 
3AMD720NUIHe 
QMD2000 365 
aMD?200U-464 
RMDPNNI 3K? 
aMD20N03468 
QAMD2NH0 SHO 
BMD? 0NU3T0 
amMOeNNU3T) 
amMD2900372 
9m02900373 
AMJANNUIT4 
gu 2nn0 37S 
amdenn03746 
QMDaDN0377 
AMD2INU3S7TS 
amOrn0379 
aM9PN009RN 
aMDADNU IR) 
am02)00382 
QmMD?P0NUIAY 
auNron0sHe 
BMD anv IAG 
aMDednl 386 
BMD2DNUSA? 
3mM0700034H 
aman IRD 
amMDd2e)n039n 
amMD2000 391 


_8MDANN0392 


aMD20N00394 
RMDAIHU 304 
AMDANAUIVE 
AMOADNVSGA 
amMD?290U 397 


a 


ete ata ee ape ee 


080276 


333 RETURN 


SJBPROGRAM LENGTH 


000112 


92 


Ta CASE SELECTION 


000277 1000 FORMAT ()3a60Ac) 
600277 2008 FORMAT(SINIILCEGAL OPERATOR O% HELAILON 9 Aee SH 
In CARD: PROGHAM SKYPPED TO NEAT CRUDE EM.) 

069277 END 

SUBPRO SHAM LENGTH 

00061f 

coo. 

FUNCTION ASSISNMENTS 

STATERENT ASSTONMENTS 

1 = 990115 2 = vo0Tri 3 = 7100125 ry 

5 - 000135 6 ~ y0ljal li = u0O1)¢ le 
“43 * 900227 14 = vuv0y33 15 = 000137 is 
20 = 900)%4 ee = uvdnre al = «0gn3? e4 
29 = 900067 50 = y00I167 > = 900031 ing 
19) = 000166 199 = yO017» eva = vbvo2n3 2? 
223 - oovl/7 361 = y0U2)7 302 = vsvo230 3046 
305 - 60026) 306 = vdlens 307 - b0ens 3a 
311 = 000214 312 = yu076) 3143 - v00222 32n 
32) => 000210 333 = u0urTte 5u0 - uon20e 50] 
i$00 » 900303 20u0 = unbaro 2¢eQ9 - wvo}03 e3eu 
BLOCK NAMES ONO LENGTHS 

= ye3elT 

VARIABLE ASSIGNMENTS 

A ~- 00606 Ale3 = gv0on3 6 - v0nends o]?4 
ce = 000603 oFe3 = yeonns I = v0060u IDdata 
rN = 909333 Is = yv0seni fy) = vv0e0e2 K 

KK = 0N06U7 KA) = vhdane NTAPE = vvONN6 NIcst 
SKY = Q0UNUVOCOL] WK = yoone) x = u00332e 

START OF CONSTANTS 

0060301 

START OF TEMPONARIES 

000322 

START OF INDLRECTS 

000330 

UNUSED COMPILEN SPaceE 

021300 

SUBKOUTINE CONVICC. SY¥M_ KA) 
CCONV SUBROUTINE CONV FOR BMVQ2ZD (3000 FORTRAN VEOSTOa) 

000006 DIMENSION CL(19005) 6 SY™(65) 9 Criak16) 

000006 TYPE INTEGER Ce Sys Cran 

000006 v0 2 é B |e KK 

000007 NECODE (6e100CC(1)) (CHad(J) Je} 90) 

000027 LO FORMAT (6R1) 

000027 00 1 J ® ly & 

000033 CHAK (J) BSYM(b5eCHAQ(Y)) 

080036 1 CONTINUE 

000040 2 ENCODE (oOe2MeLC(1)) (CHa (J) euB)eod) 

0900066 20 FORMAT (6A}) 

000065 ~ RETURN 

000066 ENO 


090132 
000323 
000343 
ovdnse 
00048} 
Qv0}73 
ougrr2 
000270 
vuvel? 
0001460 
0U00el 


000004 
O42eHILVi 
v00604 
000900 


aMO20N0 398 
RM02D90399 
muo2p0dson 
amMdz0n0e0) 
aMD?P0NVenr 


auDrnnben? 
aMDann0ens 
am Annvens 
aMD20N040K 
3M02000607 
aMOdaNn denn 
aM AnNGen9 
BMDeNROaTA 
amoanndsy) 
RMD AHHOa LT? 
amMann0412 
RMDANNNS Ls. 
AMD? 900615 
B8MD?20N0614 


oc 


COnv 
FUNCTION aSSISNMENTS 


STATEMENT ASSIGNMENTS 


i 


Fe OR Oe RR A «ST te ten th, re Ae ne ee Mt mn 


93 


| 40 = 00uvel) 20 = v00073 
j BLOCK NAMES and LENGTAS 
VARIABLE ASSIGYMENTS 
CHAR = doUulU2 1 = yoo J - nuvotil 
START OF CONSTANTS : 
000070 7 
| START OF TEMPOXARIES . gonucist 
000075 yot ® 
j 
i START OF INDLRECTS 
| 060106 
| UNUSED COMPILE SPace 
| 022600 
SUBROUTINE HATTY (agt) amM070n0417 
cPalfty SUDKUOUTLINE PATTY FOR devdeu MECIMPT  atTian VATE Hecb=4% RMOADNVETA 
| Cc : RMDANAVELY 
i 000005 DIMENSION A(135¢135) oNN(A) amyase0ean 
| c amM07900462} 
000005 IT#1 AMD2DNUee? 
000006 KK BY AMD20N062 3 
000007 w1BIT RMOPINV4G 24 
000010 Ken MIN (Bei) aM0200U405 
000013 5 AKsAKes BMDANN0e 26 
o0v0Ls IF (N=KK) 39 39% BMD ANNUSAT 
000017 & TTsl[Tel AMDPDAVGAR 
J 900021 GO To 5 RMD?INUe29 
; ‘ 00002) 3.00 $0 JaAsielT AMD 20004 
| 000023 LLLEK2en1¢1 aM020N04631 
000026 LLau AMDZ9NUS 32 
000027 )0 40 JuaKlened QAMI7NN0433 
} 000030 LLALLel BMD2NNUS Te 
| 000032 40 NN(LL) BIJ 3M02900425 
| 000036 WRITE ( 69300) (NVC TT) odd Fl OLE? BMD 2000424 
000050 DO 10 T®loN ; BMDAINVG I? 
' 000053 10 WRITES Heh Te(A(Loy) even lene) am0900438a 
i 000076 Ki ak2el 3M02000439 
i 000100 K2an 167 RMD70004460 
| 000101 K2s MENU (Keen) RMD71NDG4) 
: 0001046 20 FORMATCIA Lasky] o4e7Fla.e) ; QMD200066? 
; 000104 300 FORMAT CLHOBASSHCUL, 7(LOXPSNCOL 2) SOAS T Jel ( VhAst a) oe shm RUWS/) Onn 04s3 
00vlL0s bo CONTINUE "512900466 
| 0G0lC7 RETURN =: pN00645 
000107 eNO ar! 2DAVebA 
SUBPROGRAM LENOTH . 
000157 


“ ‘ 
Aan nrwememerre i; os = CRE 
PATTY 
FUNCTION ASSIGYMENTS 
STATEMENT 45S) SNMENTS 
3 ~- 9000¢) 4 =~ g00n17 be) ° vdols 2n » ov00113 
300 . = 000317 
BLOCK NAMES aNO LENGTHS 
VARIABLE ASSIGYMENTS 
1 = 000155 Bi >= y0u1S4 ii - w0Mee J = 000158 
dd = 06153 JK = yG60)50 KK = C0N a5 Kj = 0002446 
Ke 2 900)47 Lb = yo0TSe © LeL - obG15s NN ~ 000136 
START OF CONSTANTS 
000112 
START OF TEMPORAKIES 
000126 
START OF INDIRECTS paonuclBtt 
900133° not RE 
UNUSED COMPILEX SPACE 
922300 
SUSKOUTINE PLUT (NV st USNDL ee AMLy) QAMOPNNUAET 
cPLuT SUBROUTINE Fok omMp)eu NUV 3, 1966 QuQ209044H 
; ¢ 3AMD0?000469 
¢00010 TYPE INTEGEH By eAS AZOPy ADAG SO PAH ADeALO ALT a UMenr eCael je Si aM)20004654 
1 oleCCotICeoP ed HUF gh D0C QM02900451 
000010 . TYPE INTEGER CumeCcoce 3M029N0462 
000010 DIMENSION CIC (295) eKUF (4) am)ann0asa 
000010 DIMENSION. COCE INS) 6 LC (1300S) oC S10) oCM (6) 9 EK (159) oe TV(TSO) OXC LDU) ot 8MD7)00456 
£ (255) 94m (20) PAMIN (ISU) (L203 OT ELL SDD (ROT01 5) 8e YM (aR) OIC IB) QM)? 200455 
Cc a aM0?9004556 
00001C COMMON COoXeCoTKelvyeMMe Tr eSUMeC Jeb : 2M02)00457 
000910 EQUIVALENCE (CooTCofd) 6 (CS uC) AMDPNNDAGRA 
000010 - EQUIVALENCE (CEC oC) amne70ueso 
000010 NV LaNVel QM071)N0469 
oov0r2 Als(+64//9$%) 3M07)0044H)] 
000013 A2a(+6NSSSSO*) QM02000482 
000015 “A3s (+ONe eee >) ; QMDPDNNGKY 
CCCCOCCCCCCCCCCCCCCC Crane FOR Nitec TARCH 2968 COCCCUCCeccCc BMDPNNVGHG 
: Cc A4E(HhHoerern) 3M0720004645 
000016 ASE (HHe tet~) RMDANNUaKA 
000020 PL2(*6NPLUTSL) amDP0AN0eA7 
00002) A5 abu sene7 AMD 2 UN0a6A 
000023 AGSEHYAWVUT AMd?20N0469 
000024 ATROHSKYPON RMIPDNN4TO 
000026 ABBOHMLK UI AMD?PNNUSTI 
000027 AQ ZCHGFEDCH aM02000472 
00N03) 410*6HAY8765 aMD02900673 
000032 All#6H932) BMD A)NUSTA 
C BAVTELLF COMMENT INCREASE OCTAL woRD SIZE 10 au JrotTe §M)7 0004746 
000034 CMMB3I7TO0N00NUN0NYOINGULOUN 9MN20N0474 
000035 COCHS0OVOVNOULUOUCUNNLUONK AMD2NNNATT 
0009037 DM 800 JaleNV 9M029004 7h 
0000460 CALL SCALE CAMIN( LT) 0601) 100. eI PAMIN(T) oR CT) erin) aMD?P0NUET9 
000082 800 RID BRL -aMING) AmD?000484 
C GBalTELLE COMMENT INCREASE OClAL WOR SIZE TO 20 JrbITe auhanndent 


090062 
000063 
900065 
000066 
00u070 
000071 
000073 
000074 
600076 
000077 
000105 
000106 
000161 
000160 
000166 
000166 
000264 
000166 
000167 
000174 
000205 
000205 
000232 
000232 
000241 
000243 
000257 
000257 
000263 
000265 
900271 
000276 
000276 
00030) 
000305 
000312 
000316 
000315 
000317 


000323 
000324 


000326 
000330 
000333 
060337 
000342 
000343 
0003467 
000350 
000352 
000354 
000356 
000360 
0003605 
000366 
000377 
000401 
000403 
060405 
6006)5 
000417 
i 900427 
000432 
i 000440 


PRR termine meer 


31 


7/0 


ane 


70 
Lo 


CM(1) e7TOOOCVOOULUNDUVONNLR 

CM(2) 2007 700U00U00HNUENHOUG 

CM(3) R00U0T7UUQUOLOQNULOLOOS 

CmM(&) 28000007 /0u0CNOOUGDUGS 

CM(5) BONQ0CUUN7T fULOOULOROUS 

CM(6) 200000UU00U7 (OC0000000 

Cl eOHAAAAAA 

CEB oHHeoseee) 

OO 31 Isls6 

CUT) 8C le ANU eCM(7) 

[784 

ENCODE (665400 sQUF) A} PAZ RIPASPAS SAO AT NU GAG Al re AL} 

DE CODE (65440 Le RUF) CSYM(T) 9 dm] 965) 
Ka (Kel) /6e) 

FORMAT (11 A6) 

FORMAT (72A)) 

REMIND IT 

Kx 

00 $70 LI] #)] 9 PI. 

RE AD(595G00) (x01) ottl en) 

FORMAT (BALY) 

DECODE (72e2eA) Zell oie (MmM(I) e181, cu?) 

FONMAT (Abe T 31202672) 

IF (LL@CNVIMLL)) de 907! 

IF (Z-PL) 39593 

WRITE ( BolU) Thy (X41) 915d 522) 

FORMAT (eq HOEKROR UN SELECTION CatUL39G5Xe1246) 

6G TO 9790 

09 6 121926 

IF CONVomM (TD em) 3 eT e7 

IF (MM(T)) 77 Ue TU—b 

KeK*} 

IX(K) EMM (T) 

LY (CK) SLE 

IF C(ONPL*I1) @ (Re60) 870 .8009600 

KNaK 

Lisu 

LOxL1¢) 

LIS MIND (Llel5eKN) 


Keg 

00 500 Lalovll 

AmKe] 

IM(K) @EACL) 

LY (K) @£¥(L) yo 
KK ene 7#K ) 
NO 32 ValekK ank® 
IC(N) 80 gi 

KLB255 .\\ 

00 13 JJelenu 

nO 2 Jaleny 


KL&Kio)] 

IF (NL #255) 1201 rel 
READ ( ThyCye 
KL=) 


R( J) RCIC (KL) 

DO 13 Jalek 

Lalx(J) 

malY (J) 

Lea eSe(KX(LIMAMEN( LD /Q(L) MEUe 
{F (52°) 13943670 

ma (X(M) @AMIN( A) )/K IM) #100641 65S 
LF ((192—0m) Om) 19013572 

Mle M00 (Mels6) 0] 

M28 (M05) /6 


AMY2900682. 
3Mn2)006R3 
AMDADNOGRA 
BMDADHUGAS 
840200046846 
aMD20N0487 
AMDZNNOaAR 
AMD 7/)NUSRY: 
3éOeNN0eOn 
9M070N0691: 
BMD?0N0a92 
AMD250069 3, 
amD?000694 
9M02900495 
9M070004904 
3mMD2900807 
aM enndegu 
QMD02000694 
aMD?20n0509 
3mMD700U5N) 
3m07000502 
3™02000503 
Smpe0NUSN4 
3MO?P 000505 
3MD2N0US5N6 
aMO?PD0USN7 
3M02900509q 
IMD2DNNSND 
aMD2D00510 
3m07i00511 
QaMDznn0S12 
amo7pn0513 
9mD2000514 
3M02000515 
BLD20N05S14 
3MD2NNOS17 
96D 2000519 
3MD20005)9 
amD2nn0SA4 
9mM02000521 
amdenandser 
AMOPIAUSPA 
AMDAONUSA4 
QMOANNUSAS 
AMZ INVS2K 
3MO02PINUSD? 
QmMD2y00528 
amM)20NN529 
aM02000534 
349290053) 
9M02900532 
aM02000533 
amM02900534 
4M02000535 
aMD?20N05S AK 
3mM02000837 
aM07000534 
amMdAnNUS IAG 
ame 0054an 
RMDANNUSE) 
AMD? 000567 
aM02)00563 
aMD?27NUS66 
aM02000565 


0004464 
000652 
000455 
000460 
009461 
000464 
00047) 
00067) 
000476 
600477 
000506 
000533 
000515 
600520 
900522 
000526 
060533 
000534 
00054) 
000542 
000546 
000546 
000550 
000552 
000553 
900555 
000560 
0005646 
060566 
000567 
000572 
000574 


000702 


000702 
000704 
000705 
000710 
000711 


001137 
PLULT 


2 

8 
13 
32 
60 
70 
301 
510 
701 


40 


4} 


51 


52 


5n 


56 


590 
S70 
900 


3914519 (N2¢ 172-18) 

IF (M1—1) 60960.61 
MUGECM4e AND e eNOT ICC (43) 

LF (MUG) 6196306) 

IF (CC{ M43) ) 13064066 

cc (m3) SCC (MS) sANDe eNOT eC MOURECCE 

GO TO 13 

MUGBCM (MD) eANU. eNOT.CE (M3) 

1F (MUG) 16913916 

IC(M3)MISIGN (TAMS (10643) )egC(ms) 9 £0 (M3)) 
CONTINUE 

RE@INO [T 

CALL CONVICCeSYMoKK) 

paloltHe) 

VO 40 Nal,101 
C(N) aP 

OO 41 Nebo lls 
C(NI 80% 

vo 50 Jelek 
Lair) 

mely (J) 
QeAMIN(M) 
BRIM) A106 

vO 5S] Nalell 
6(NI sq 

QzQen 
GwAMIN(L) RIL) 
DERILI 7506 

OO Se Nsle5l 
KCN) @Q) 

QeQ-) 

WRITE ( 


Lalo lOdd o (XCRU0) oC (KN0) o LDO (Ned) oNBKRU09A67 953) 9CIKON, 


S51) (CUM) ONBL OTOL) a CTIND oNB@L ebb eC) 9 (TIN) oNBA 040.) 
FoemMaT (LIA 
L/ALBAGLOIAL SDL ISIS Ve Be LAPALS IAAP AD,Al 
CoS5F2003/ TKO 2N.3) 

[LF (LI=KN) SlusdAus SA. 

n20 

CONTINUE 

RETURN 

END 


SUBPROGRAM LENGTH 


FUNCTION ASSIGYMENTS 


STATEMENT ASS] ONMENTS 


= Qnu7lo 3 = 000243 b = 000263 ’ 

= o00u2%6 9 = vuerarw lv = v0OTT , 

«~ 000506 14 = C00 FK0 10 = vd0er? 35 

= 00344 40 - v00827 cy! - 006635 56 

~~ 000655 41 - v0v6T7) 63 = vbesal 66 

= v0usI7 71 = 0432 vf - v002e1 Yo 

= 000352 302 = ooOuaTT aud = 00741 41 
= 09000315 Tu = 0070S 540 - vloTos on 
= 000813 Tlo = 6009S yun - vodoTley Sonu 


= 0327 


BLOCK MAMES AND LENGTHS 


2 


ho DH) 9L 9 (TIN) ONE] L102) 9 CT UND NBO, 1092) og (COND GN 


KU0BLe 


VARI AAL EAA SHAVORIASLEL3/TT/2K eo FILS, SoS Ae A/T NSP EU, 3 
)/LIM, LOLAL/AXgh 15,2 


Quo77) 
000346 
geo 00 
Ovlole 
000464 
VU0164 
0007463 
000312 
0U0766 


BMO?0N0S46 
3M02000547 


gm07000546R° 


3mM02900869 
340720005506 
3m02900551 
3M02Nn0552 
4m0?0N0553 
3mD2000556 
am02000585 
9M02000554 
am02900987 
3¢02000558 
3m02000589 
3MD2N0U56N 
ax02000561 
amO02000S62 
aM02000563 
3mM0700US#4 
9607900565 
9M0200054646 
9m02000567 
8402N005HA 
amM020n05469 
AMDPDOUSTO 
340290057) 
RM02000872 
9M02N00573 
RMDANNUS 74 
amMOenn0S76 
AMDADOUSTA 
3M02000577 
aMD2D00578 
aup?700874 
3M020N05RN 
9M0290005R) 
QuDannysSRa 
8402000983 
AMDANDUSAS 
AMDADOOSAS 
QMO> INOSARH 
QM020N0587 


cp GT, eS NAS uses PSST Pg NN 


Fn a nt ar er a et eR AE A A Nn eR -nRR S 


VARIABLE 
Al 
a3 
ar 
c 
CJ 
ce 
I 
1x. 
JJ 
KL 
LL 
mM 
m3 
PL 
x 


START OF 
000713 


START OF 
001032 


START OF 
001046 


ASSIGNMENTS 

* 901053 aly - GOlLNES Ali = 001065 a? = 0020866 
001055 as = wyoln7e ao - ooLos? As = 091060 
00106) ae = udlnae a? - 01063 oe o PO1NT3 
031543C01 CC = v0GNNCLO1 CLC © v9l1O> Cte e 031643001 
032756C0i CM = 32746600) CMM - 092106 C1 e 0014066 
001057 iB) = o0y725 WT) = vd0N00CO! Ha ~ O01}22 
001107 rc = oodnouCcRt TA - vublli« 1T -~ ovi312 
03214¢2C01 IY = o32370C0) J = vdvli2ze Jo = 04275601 
vo11es JJJJ = 001110 K - w01}13 K« ~- 004323 
volles KN = 0011317 KUO © v011360 L - 001122 
vollis LO = vfllal Li - vv01}2u M - 0011327 
V32C16COl MUG = 001733 ml = 901130 Mm? = 901)331 
001132 N = vdl)lb NV} = 01196 Pp = 001072 
001056 Q = v01}3 stm = 32659Cn) | = 3264241 
O3L3IE “01 Z = v0in70 

CONSTANTS 

t 
TEMPORARIES wo 
Q 
ar® 
INDLRECTS Qo 


UNUSED COMPILEX SPACE 


017500 


000017 
od0017 


000017 


000017 
000036 
000040 
000041 
000043 
000046 
000051 
000054 
000057 
000062 
000064 
000067 
00007) 


0002135 
000117 
000121 
000122 
000132 
000133 
000135 
000140 
000140 
000161 


Cc 
Cc 
c 
Cc 


305 
310 
315 

5 
4 
§90 


lo 
1 


8 


20 
M1 


le 


SUSROUTINE TRNGEN (DOTA TN eNVOeNVUR! Oy TSAMDe CACC yNNEWMeLLCUDE? 


1} LEV Ag BSNEWe MERHY 9 N) 
TRANGEN 


SUBROUTINE TENGEN POR BpMuveU VUVe MFO 135 lore 


UITMENSION DATA( 150) sNNEwWA (150) SLECUUE (150) ot LYM 1150) FRONEM (TSU) 
DIMENSION VECTOR (1949134) 

COMMON VECTUR 

ASNF (AA) BATA CKAX/SORT (LeUeKxAee)? 

1 TEMaNn 

SAMPaNUUOTA 

00 3 JBieNVYs 

NE WAMNNEWA (J) 

LCVOE@LLCUDE (J) 

LVaaLLVAl(y) 

HNEwaBSNEw (J) 

1F(LCODES10) &,a05 

NE qa NEw 

p2vaTal(tva) 

IF (LCOVE #41) 900017023 

GO TO C109 2093 BU SD eb e MUG HU WVosUD eT LU ve LAV l Welaty 


11900160) eLCvvet 


IF (U)99—79A 

VATA (NEWA) U0 

60 TO 3 
DATA(NE WA) saWURT (LU) 
gO TT) 3 

LF (U0) 9Vellele 
DATA(NEWA) 1. 
60 10 3 

VATA (NEWA) BSURT 
60 To 3 


(VU) +SURT (Veo) eQ) 


aMD200058H, 


amMD2000588 
RMDzNN0HaN 
am02ijnusa) 
aMD20nUSY?P 
am02nn0saa 
3m0270005946 
BMDZNNNSIS 
3M07 INDIA 
au02)00597 
aM920N0594H 
amD29n0590 
amopo00600n 
am02nnven} 
amDenn0on? 
9MD20N06N3 
RMQ070000%4 
guDandnuens 
aMQenNn0b0n 
amOann0oor 
aMN2A NN QANN 
aAMOANNVOEDYS 
amMDannveln 
ampajndoll 
SmuOrNnvoI? 
amM02790061 3 
aMD20n0014 
amDann0e1s 
aMOzINVOIA 
aM029NU617 
aAmMip2NNGe 1s 
au02nnug)a 
anDz000020 
am0290062) 


000161 
000163 
600174 
000174 
000204 
060205 
090207 
000212 
000215 

- 000215 
000220 
000230 
000230 
000233 
000237 
000240 
000242 
0900258 
000255 
000257 
000272 
000272 
000275 
000277 
030313 
000313 
000314 
000317 
000320 
000323 
000326 

| 000327 
000330 
000331 
000336 
000337 
000344 
000344 
00035) 
09035) 

} 000356 
000356 

| 000360 
| 000365 
i 000365 
000367 
000372 
000376 
900377 

| 000403 
i 000400 
00041) 
9004)) 

i 000012 
H 000415 
i 000417 
: 600025 
| 000433 
H 000437 
000486 
| 000445 


aU) 


“2 
1604 


IF (0) 99999_h4 

DATA (NEWA) BAL0G (1D) #,4362904819 
60 Tn 3 

DATA(NEWA) SEXP (D; 

gO To 3 

IF (U)99079)7 

TF (V91.60) LB919.99 

DATA (NEWA) 23616159265/2,0 
69 To 3 

AwSORT (0) 
VATa(NEWa) mASNE (a) 

60 To 3 

AaDds (SAMPol1 60) 

RAO] 20/7 (SAMPO) 0) ‘ 
IF (A)99,23,e4 

LF (69907 s27 
DATAUNEWA) E4ONE (SORT (8)) 
60 To 3 

1F (4) 99e26e0e4 
DATA(NE@A) BASNE (SQRT (A)) 
60 TO 3 

awSeRT (Aa) 

BeSQRT (8) 

DATA (NEWA) BASNE (4) ASIF (35) 
60 TO 3 

IF (D) 31999314 
DATA(NEWA) al eUsi 

GO 10 3 

DATA (NEWA) ei eonew 

60 TO 3 
DATA(NEW2) x #bNEW 

G0 To 3 

LF (VU) 3307933 

DATA (NEWA) sU**#BANE w 

GO To 3 

DATA (NEWA) BUSDATA (NE we) 
GO T0 3 

DATA (NEWA) BU—DATA (NE wt) 
GO fo 3 

DATA(NEWA) sOMUATAINE wh) 
GO To 3 

IF (DATA(NEWO) ) 34999, 34 
VATA (NEWA) SD/0A TA (MF Wb) 
GO TO 3 

HNEwanEwy 

LF (Des New) Tobie il 

IF (UeDATA(NEWB)) Tellell 
IF (vu) 3950303 

LF USTGN (10eU0M)) §f493,3 
WATA (VEWA) avew 

CONTINUE 

GO TO 42 

LCasEa-999 

LE IMERRY@ J) 4020401 .402 


mE RKYSJ 

WAITE ( ReleNady 
WRITE ( hy 1405) 1TEM 
wRITE( 6.140K) 
JSAMPwiSAMP=] 

RETURN 


an02900622 
am02000023 
amD 2000606 
aM070006265 
amD2000624 
amD2000627 
amMD20n0628 
aM02009629 
9M0700063n 
3m0200063) 
9072000632 
34020006033 
3mM02N00636 
8402000638 
AMO? 000694 
AMD?2DN0627 
Bm02900638 
BM02000639 
RM02000646" 
AMO?DN06e) 
RM02000662 
RM02000663 
3MDPNN06E6 
907000645 
8MD7N00644 
amD2nn0ba7 
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IA. Tape Input 


Most of the statistical programs have provision for in- 


put of data prepared on tave by mea:.s of an Alternate Input 
Tape. Use of tapes for data provides a very ccapact storage 
of data files. The tape for one data file may be reread sev- 
eral times tc allow several analyses from the same data file 
by rewinding the Alternate Input Tape. Of course, data in- 
put may be from data cards if preferred. 


IB. Standard Data Input 


The form of Standard Data Input is given below: 


Variables 
7 x) Xo X3 oo we Ue x. 
1 
2 X01 Xo X53 ear Xo 
Cases 3 X34 X39 X33 ites X3p 
n x 


The headings Kye Xqe coer X represent variables, e.g., 
age, sex, weight, etc. Each row in the table represents a set. 
of corresponding values of these variables, e.g., the age, sex, 
weight, etc. of a given individual. The entries Xi in the 
table are called data values, the whole array of these numeric 
values is called the data matrix, each row of the data matrix 
is called a case, and each column is called a variable. 


The Standard Data Input is keypunched case-wise. That is, 
all the data values of the first case are keypunched in order 
on one or more cards. Then starting on a new card the second 
case is punched, etc. Each case must have the same format 
(see Section II-B). This means that from one case to the 
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II. 


ITA, 


next, each variable value must occupy the same physical lo- 
cation on the card into which it is punched. 


In general, not all of the fields on a card will be con- 
sidered as part of the data matrix. For example, identifica- 
tion fields such as the subject number are seldom included. 
The desired fields are selected by means of the Variable For- 
mat Card specification (see Section II-B). 


Preparation of Program Control Cards 

The statistical programs are writcen in a general form so 
that a wide variety of problems combined with optional compu- 
tations may be handled by each program. The user specifies cer- 
tain parameter values, optional computations and opticnal 
output, the form of the data input, etc. on Program Control 
Cards. Standard program control cards which are used in many 
programs are described in this section. Instructions for the 
preparation of other control cards specifically required for 
an individual program appear in the program description. 


Unless otherwise stated, each numeric field of a control 
card should be punched without a decimal point. The decimal 
point is assumed to be at the extreme right of the field, and 
blank columns will be interpreted as zeros. Thus, for example, 
a field including Columns 1-6 (1 through 6) which contains only 
a 1 in Column 5 will be interpreted to mean 10. In general, on 
control cards numbers should be punched in the rightmost col- 
umns of the field. This is called right-justification. 


Transgeneration Cards 

The term transgeneration is used to include transforma- 
tions of input variables and creation of new variables prior 
to the normal computations performed by the various programs. 


The transformations described below are performed on the 


values of the variables in each case. In these examples, the 


symbol x, will denote the ith variable as weli as its value. 
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Examples: 
10g, X4 * X4 log, 9X4 replaces X4 
c c 
Xe 7X) X. replaces X, 
X, + x, * X5 Xy + X3 replaces Xo 


By successive transformations, more complicated rela- 
tionships may be obtained. For example: 


(1) To replace xX, by x,“ + x? four transformations are 
required: 
Variables as they are stored at each step 

Transformation 

2 + 
ey ay 

2 
Xy Bc X3 
Rey ig ae 

vk wits 


In this example, it can be seen that the original values 
of X, are irrelevant. Actually the variable X, may be a dummy 
variable introduced by the program specifically to provide 
capacity for creating new variables by transgeneration. Dummy 
variables may be required for intermediate storage in order to 
effect some transformations. 


3A 


(ii) 


To replace xX) by exp (-1/2 x17) three transformations are 
required: 
Transformation 
2 
“1/2 X, +X -1/2 x,” 
1 “1 1 
exp (X,;) +X, exp (-i/2 x,7) 


(iii) To replace X, by X, + 109, 4 X4-X, +100) four transforma- 


tions are required: 
Transformation 


Xy-X, 0 XQ 


lag, 4 (X4-X, +100) 


The transformations are performed in the order in 
which the Transgeneration Cards appear, so that, for example, 


the two transgenerations 2X, + xX, followed by Xj,-2 + Xy wiil 


+ X, will 


result in 2X,-2, whereas X 1 1 


1 
result in 2(X,-2). 


172 + X) followed by 2X 


TRANSGENERATION LIST 


Notation to be used in the following transgeneration list: 
i, j, k are variable indices (need not be different) 
c is a constant 
a. a 


3° ay eee are constants 
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n is the number of cases, or sample size 
3 
SA 


The standard deviation s.= 


Code Transgeneration Restriction 


01 X,* x, X, > 0 

| 02 VK; + /x,41> x, X,; 2 0 

| 03 log, 5X, * X, X, > 0 

| 04 = *X, = 

| 05 arcsin /X, + x, | 0<X, <1 

06 arcsin VX, 7 (atly ‘ arcsin /(X;+1)7(n+1) + x, 0 < (Xm) <1 
: 07 1/K; + X, x, AO 

| 08 =X; tc* X am 

| 09 X;c * x, - 

| 10 x,° +x, 20 


11 a Ba ie Xx. a 


12 Xi-X, oo X. 


13 XX, sale oh - 
14 X,/X, 7% x #0 
15 Tf X, > co 1 +R; - 


otherwise 0 * xX, 
16 Xf a, 2 Xye i 7X - 


otherwise 0 * xX, 


i? icg,X; * x, x; > 0 
F 

p vw + - 

; 18 xy x x, 

d 6A 


Code Transgeneration Restriction 


19 X,/8, + x, i 
21 cos x; * x, | - 
22 arctan XxX; * x, - 
23x. *5 + X.> 0 
i x, i 
Xs 
24 ¢ * X, c> Q 
25 x,* X, = 
26 c +X, (Leave code i 
blank) 


27-39 Not defined 


40 If X; * a) OF A) OF a, .+-, ay, then c + Xi 
otherwise X, remains unchanged. 

41 If X, is blank, then c + XF (X, # -0)* 
otherwise x, remains unchanged. 


*Note that in reading numeric fields, a blank 
field and -0O are equivalent. 


42 If X; = @, Or a, OF a, «++, ay, then x, > X3 
otherwise x, remains unchanged. 

43 If X, is blank, then x, > 7 (x; # -0) 
otherwise x, remains unchanged. 


When a violation of a restriction in the right-hand column occurs 
during transgeneration, the program will print a diagnostic mes- 
sage. Most programs will proceed to the next problem, if any. 


TA 


aes eg ae ae te 
L100 = ARENSON TP TIES OPS SEE 


vi dt dts Cee o eeaclaeea eee 


ete ae eng ee A EFI AAEM STE OIE 


Some programs will delete the case where the violation occurred 
and continue the computation. Other programs will screen all 
the input data from additional restriction violations before 
proceeding to the next problem, if any. 
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1. Standard Transgeneration Cards 


Standard Transgeneration Cards are used with programs which 
use Standard Data Input (see section IB). Let p denote the 
number of variables in the data matrix, m the maximum num- 
ber of variables allowed by the program for any problem and 

q the number of variables added through transgeneration. 

Any of the variables po coer XH, may be used in transgener- 
ation. The initial values of the first p variables are read 
from the input data file (Data Cards or Alternate Input Tape). 


The initial values of the remaining m-p variables are left 
over from previous caiculations. After transgeneration 
action of a particular case, the values of the first ptq 


variables for that case are used as the values of the trans- 
generated variables. If the ptq variables required for the 
computation are not the first p+q, they must be relocated. 
This may be done by using transgeneration code number 25. 
The numbers p and q (q may be positive, negative, or zero) 
are specified on the Problem Card. The indices i, j, and 
k from the transgeneration list may exceed p or ptq but must 


never exceed m. 


Card Preparation 
Col. 1-6 TRNGEN (Mandatory )} 


Col. 7-9 Variable index k 


Col. 10,11 Code from transgeneration list (restricted 
by availability in particular program) 


Col. 12-14 variable index i 

Col. 15-20 variable index j or constant c 

Col. 21-25 Blank 

Col. 26 Number of a,'s for transformation 40 or 42 
Col. 27-32 a) value 

Col. 33-38 as value 
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IIB. 


ase 


Col. 63-68 ay value 


The constants c, Aye ose, a, are punched with a decimal 
point if used with variables which have an F-type formt 
and without a decimal point if used with variables which 
have an I-type format (see Section II-B), 


The Standard Transgeneration Cards for the three samples on 
pages 3 and 4 are: 
(i) TRNGENOO01100012.0000 

TRNGENO03100032.0000 


TRNGENO0511001000003 
TRNGENO0501005000000 


(ii) TRNGENOO1100012.0000 
TRNGENO0109001-0.500 
TRNGENOO01 04001000000 


(iii) TRNGENO0412004000003 
TRNGEN00498004160.00 
TRNGEN00403004000000 
TRNGEN00411004000002 


Variable Format Cards (for Input) 

The word "format" usually refers to the arrangement of 
information keypunched on a card. The format of a data card 
is a sequence of fields (variables), each of which occupies 
one or more columns. For the computer programs a format is 
a set of specifications according to which information is 
read into the program from punched cards. The specifications 
tell the program which parts (or columns) of the card to skip, 
which parts to regard as all one number, and which parts to 
regard as several numbers in a row. For instance, it is the 
format which tells the program whether a card punched 
"345890021" is to be read in as "34.5, 890.0,.021", or "34., 
9002.1", or "589.002", or "3, 4, 58, 90, 0, 2, 1", etc. It 
does this by givi-q the program a sequence of specifications 
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which indicate the size of a field and the method of han- 
dling that field (i.e., skipping ic, entering it into the 
computer as a whole number, entering it into the computer as 
a number with two decimal digits, etc.). 


The format also tells the program how to read in a cer- 
tain set of cards when more than one are required to contain 
all the data for one case. In most programs the format de- 
scribes the variables for each case. Each successive case 
is assumed to have the same format. 


If the formats for the programs were fixed in advance, 
all data would have to be punched on cards in the same way 
for every study. Since this is not usually convenient, the 
statistical programs have been written so that the user may 
vary the formats according to his preference for a particu- 
lar study. For this reason, they are referred to as "vari- 
able formats". The program is informed of the format which 
is to be used by Variable Format Cards. The user must spec- 
ify on the Problem Card the number of cards used to keypunch 
the variable format. 


In addition to providing an economical method of pre- 
paring data input cards (by defining fields to be as small as 
possible, or "packing" the data), the variable format permits 
considerable freedom in controlling data input. For instance: 


. It allows the user to select for each case only 
those cards which have fields of interest. (See 
Examples ii, iii, vi helow.) 


. It allows the user to select only those fields of 
interest from among the fields of each card. (See 
Examples iv, v, vi.) 


. It allows the user to scale the data input, i.e., 
shift the decimal point. 
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| A complete description of formats can be found in 
| FORTRAN programming manuals such ae those available from 


IBM representatives. The features commonly required for 
the statistical programs are described below. 
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1. F-type Variable Format 


The F-type format is the most frequently used in the 

statistical programs. It is required when the deci- 

mal point is keypunched on the card or when the dec- 

imal point is to be placed by the program. All data 
input values must be signed (+) or unsigned numbers 

with or without a decimal point punched. | 

: 

| 

| 

| 


Specifications: 


(a) “nFw.d" F is the floating point (decimal) 
indicator; n is the number of fieids of 
width w (includes sign and decimal point if 
punched); and d is the number of digits to 
the right of the decimal point if the deci- 
mal is not punched (0 < d <w). If the 

decimal is punched, d is ignored. If nv is 
not specified, it is assumed to be l. 


(b) "sPnFw.d" P is the scale indicator: s isa 


scale factor (explained below); and n, w 
and d are defined in specification (a). 


| (c) “mX" (alphabetic X) xX is the skip indi- 
' cator, and m is the number of columns to be Ce 
skipped. Nae 


| 
| 


{d) “/" (slash) indicates "go to the next data 
card." 


Depending on its location in the format state- 
ment, the "/" will either direct the program 
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to go immediately to the next catd (ignoring 
any further information on the current card) 
or skip one card altogether. For example, 

if a format begins with "/", the program wili 
automatically skip the first card, read the 
second, skip the third card, etc. Ifa 
format ends with "/", the program will auto- 
matically read the first card, skip tie sec- 
ond, read the third card, etc. (See Examples 
ii, did), vy vis) 


"//" indicates "go to the card after next." 


Two slashes "//" will direct the program to 
skip two cards. Any number of slashes may 
be used. 


The format is keypunched beginning with a 
left parenthesis, a sequence of specifica- 
tions, and closed by a right parenthesis. 
Specifications (a), (b), and (c) are followed 
>y a comma, except preceding a slash or 
right-hand parenthesis. Blank columns within 
the format aze ignored. Columns 1-80 may be 
used unless otherwise specified. 
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(4) 


(ii) 


(iii) 


(iv) 


(12F3.0,F4.0,11F2.0), punched in the first 20 columns 
of the Variable Format Card, will describe 12 three- 
column fields followed by 1 four~column field, followed 
by 11 two-column fields. Each data card will be read 
according to this format. 


(12F3.0,F4.0,11F2.0 /), punched in the first 22 columns, 
will describe the same :ields as in Example (i), but 
will also instruct the program to read the first card, 
skip the second, read the third card, skip the fourth, 
etc. 


Two slashes at the end 
eee, LLF2.0 //) 


will instruct the program to read the first card, skip 
the next two cards, read the fourth card, skip the next 
two, read the seventh card, etc. 


(/12F3.0,F4.0,11F2.0), punched in the first 21 columns, 
will describe the same fields as in Example (i), but 
will also instruct the program to skip the first card, 
read the second, skip the third card, read the fourth, 
etc. 


{Two slashes at the beginning 
(//12F3.0,.0. 


will instruct the program to skip the first two cards, 
read the third, skip the next two, read the sixth card, 
skip the next two, etc. 


(10X,F6.0,2X,2F3.0), punched in the first 19 columns, 
will instruct the program to direct entry of data from 
each card as fo).lows: 
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(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
; (7) 
(8) 


(9) 


Skip 10 columns. 

Pick up a six-digit field in Col. 11-16. 

Skip 2 columns. 

Pick up 2 three-digit fields in Col. 19~21 and 
22-24, 


(v) (5X, 2F6.0,F1.9,3X,F5.0 /5X,F6.0), punched in the first 
32 columns, will instruct the program to direct entry 
of data from each pair of cards as follows: 


Skip 5 columns. 

Pick up 2 six-digit fields in Col. 6-11 and 12-17. 

Pick up a one-digit field in Col. 18. 

Skip three columns (Col. 19-21). 

Pick up a five-digit field in Col. 22-26. 

Go to second data card. 

Skip 5 columns. 

Pick up a six-digit field in Col. 6-11 of second 
card. 

Repeat for each pair of data cards. 


(vi) (5X, 2F6.0,F1.0,3X,F5.0 //5X,F6.0,/F4.0,2X,F1.0), 
a punched in the first 47 columns, will instruct the pro- 
cad gram to direct entry of data from each set of four cards 


as follows: 


(6; 
: {7) 
(8) 


(9) 
(10) 


(11) 
(12) 
(13) 


(1)-(5) Sample as Example v. 


Skip second card and go to third card. 

Skip 5 columns. 

Pick up a six-digit field in Col. 6-11 of third 
card. 

Go to fourth card. 

Pick up a four-digit field in Col. 1-4 of fourth 
card. 

Skip 2 columns. 

Pick up a one=-digit field in Col. 7. 

Repeat for each set of four cards. 


SCALING 


Scaling may be indicated by using either the "nFw.d" speci- 


fication or the "“sPnFw.d" specification. This specification 


will not often be required. 


"nFw.d" Specification: 


When the decimal point is not punched, the d of the above . 
specification instructs the program to divide the whole num- 
ber picked up by 10%, For example, F6.1 will specify that 
the number picked up in a six~-column field be divided by 10, 
F6.2 will specify division by 100, and F6.6 will specify 
division by 1000000. 


Stored number = punched number /109, 


When the decimal voint is punched, d is ignored. 


Examples: 

Punched Format Stored 
Number Specification Number 
2468 F4.0 2468.0 

3691 F4.1 369.1 
4810 F4.3 4.810 
4911.32 F7.0 4911.32 
4911.32 F7.4 4911.32 
172115 F6.6 0.172115 


The entire format for these specifications might be punched as: 


(F4.0,F4.1,F4.3,F7.0,F7 -4,F6.6) 


"sPnFw.d" Specification 


Whether the decimal point is punched or not, the s of the 
above specification instructs the program to divide the num- 
ber picked up by 10°. For example, 2PF6.0 will specify di- 
vision by 100, -2PF6.0 wil! specify division by .01. 


Decimal point punched, 


Stored number = punched number /10°, -8 < 8s < 8 


Decimal point not punched, 
s+d 


Stored number = punched number /10 


Examples 


Punched Format Stored 
Number Specification Number 
7432 1PF4.0 743.2 

74.32 -5PF5.0 7432000.0 
7432 2PF4.3 -07432 
7432 1P2F2.0 7.4 & 3.2 
7432 2P2F2.1 0.074 & 0.032 
7432 OPF4.1 743.2 

72432 -SPF5.0 7432.0 


The entire format for these specifications might be punched 
as: (1PF4.0,-5PF5.0, 2PFP4.3,iP2F2.0, 2P2F2.1,0PF4.1,-3PF5.0) 


Mixed "nFw.d" and "sPnFw.d" Specifications: 


Once the sPnFw.d specification has been used, it will hold 
for all Fw.d specifications to the right of it until another 
sPnFw.d is encountered. 


If the sPnFw.d specification is not necessary to the right of 
its occurrence within the variable format statement, then 
OPFw.d should be used for the next specification to the right, 
which will hold then for the remaining specifications. 


Examples: 
(vii) (2F3.1,2X,F4.2,-6PF3.0,F4.0) 


The program interprets the last field as -6PF4.0. 


The program interprets the last field as F3.0. 


(ix) (1LPF3.0,F2.0,O0PF3.1,F4.2,F4.0) 


The program interprets the second field as 1PF2.(G. 
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I-type Variable Format 


This format is required for programs designed to process 
only integer values. The specification is “nIw", where w 
is the width of the field (includes sign if punched), and 
n is the number of fields (assumed to be 1 if not punched). 
All data must be signed (+) or unsigned integers with no 
decimal point punched. Examples corresponding to those 
given for F-type are: 


(i) (1213,14,1112) 
(141)  (1223,14,1112 /) 

seeg lTI2./7) 
(iii) (/1213,14,1112) 

(//1213,... 

(iv) (10X,16,2X, 213) 
(v) (5X, 216,11,3X%,15 /5X,16) 
(vi) (5X,216,11,3X,15 //5X,16 /14,2X,11) 


No scaling is permitted with I~type format. 


A-type Variable Format 


This format is required for programs designed to process data 
with alphabetic, numeric, or special characters, or combina- 
tions of these. ‘The specification is "nAw", where w is the 
width of the field, 1 < w < 16, and n is the number of fields 
(assumed to be 1 if not punched), Each specification of a 
field results in a computer word consisting of exactly 10 
characters. When w < 10 the characters are positioned in 


the left of the field and the remaining characters are 
filled in with blanks. The following examples illustrate 
certain rules. 


{ 
| 
Punched Format Stored* 
| Data Specification Data _ 
12.0 A4 12. 0bbbbbb 
AGE A3 . AC Zbbbbbbb 
CANCER A6 CANCERbbbb 
$ Al Sbbbbbbbbb 
X+Y=A AS X+Y=Abbbhb 
DX, @R,DY A8 DX, SR,DYbb 
AGE AG bDbbAGEbDbbb 


*b indicates blank character 


| 

| 

| 

IIC. Variable output Formats | 

| All the gtatistical programs use Variable Format Cards | 
to describe the input data; a few require their use to de- | 
scribe output data, that is, data to be princed, punched, or | 
written on tape by the computer. The function of the Variable | 
Format Card is the same for input or for output: it is a de- | 
scription of the data in the medium external to the computer. 
Input and output formats are identical except for the follow- 
ing minor differences: 
1. In F-type formats ("nFw.d" specifications), the decimal 
point is present (except when d=0) in the output medium, 
and a column must be allowed for it. 


2. in using the scale factor specification of the form 
“sPnFw.d", the external represe:utation of the number is 
10° times the internal number. Thus, if the internal 
number is -15.9357, a specification of 2PF9.1 would give 
"bb-1593.6" in the external medium (punched, printed, or 
tape output). 


3. Each "line" of an input format for cards or of an output 
format for punched cards must not exceed 80 characters in 
length. Each "line" of an input format for an alternate 
BCD input tape or of an output format for printing or for 
a BCD output tape must not exceed 136 characters in length. 


4. Position 1 of the printed line is used to control paper 


spacing and normally should be left blank (to produce 


single spacing) by using "1%" as the first specification 

of the format. (Thus, when printing, only 135 positions 

are actually available to contain information.) If double 
spacing is desired, it may be obtained by using "1HO" in- 
stead of "1X". 

Note: Care must be taken to allow sufficient width for 
the maximum size number that may be described by 
the format specification. (In describing input 
formats, this is essentially automatic because 
it is known how many columns of a card are de- 
voted to a particular number.) 

Example: Suppose it is desired to write an output format to 
print 27 signed numbers in F-type format to 4 decimal places, 
the maximum absolute value of the first 13 being less than 
1C00O and of the last 14, less than 50. The specification 
13F10.4 allows room in each number for: 1 space, sign, 3 
digits before the decimal point (< 1000), the decimal point, 
and 4 digits to the right of the point. This gives 130 char- 
acters on the first line, plus 1 for position 1 (for spacing 
control), and the first line is filled. Similarly, 14F9.4 
may be used for the last 14 numbers. A new line is started 
by the slash ("/") in the format statement. If the spacing 
control is to be a double space before printing the first line 
of the group and single spacing within the group, then the 
complete format ‘vould be (note the use of "1X" in line 2): 


(1HO, 13F10.4/1X, 14F9.4) 


If it is desired to separate the numbers by more than one 
space, the following format might be used (three lines will 
be necessary): 


(1HO, 11F11.4/1X, 2F11.4,10F10.4/1X,4F10.4) 


For further control of line spacing in output formats, n con- 
secutive slashes will produce n-l blank lines. 


IID. Finish Card 


This card will notify the program that the entire job 
is finished. The program will complete its computations 
and will return control to the system monitor. 


The preparation of this card is as follows: 


~> 


Col. 1-6 FINISH ! 


Bi AMINE “nate — a ee 


III. Preparation of System Cards 


The system control cards listed in this section are 
Gescribed in detail in the User's Manual. 


A. Job Card 


See User's Manual section 3.1.1.1 THE JOB CARD. 
B. REQUEST Card 
See User's Manual saction 3.1.1.13 THE REQUEST CARD. 


Cc. LIBCOPY Card 


See User's Manual section 5.1.1 BINARY LIBRARY (re- 


placing CSCBIN with STATBIN). | 
D. LGO Card | 
See User's Manual section 3.1.1.5 THE PROGRAM CALL 
CARD. 
| 
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a Iv. Examples of System Cards 
A. Job Card 


1. 7777,JONES,CM6G000,T109. 
CM field length = 6000C (octal) 
time limit = 100 seconds ‘!decimal) 


2 priority level = 1 (decimal) 


2. 111.1,SMITH,P5,T64,L1100. 

C fieia length = 4000 (octal) 

time limit 64 seconds (decimal) 
5 (decimal) 
1100 lines (decimal) 


on 
e 
e 


priority level 


line Limit 


oo 3. 33333-BROWN,CM50000,T85,P10,1L3306. 
— CM field length = 50000 (octal) 
, time limit = 85 eeconds (decimal) 
priority level = 10 (decimal) 
line limit 


i] 


3300 lines {approx. 50 pages) (jecimal) 


B. LIBCOPY CARD 


1. LIBCOPY(STATBIN,LGO, BMD8V) 
2. LIBCOPY(ETATBIN, LGO,WRAP) 


C. Typical Deck Set-up 


Jobeard with,CM70000, P5,T32,1L1100. 

LIBCOPY (STATBIN, LGO,WRAP) 

L6O. 

7-8-9 {end-of-record card, :.ultiple-punched in column 1) 
(program control cards and data) 


6-7-8-9 (end-of~information card, multiple-punched in 
column 1) 


woe tener mermar ness 


D. 


Source Library 


If a source listing of the program is needed, the 
source decks are contained in the COMMON file named 
STATLIS. 


The following example shows how to obtain a 
source listing of, compile and execute program BMD3R. 


Jobcard with CM77000,1L3300,T64. 

COMMON (STATLIB) 

EDITS YM (C=COMPILE , OPL=STATLIB} 

RETURN (STATLIB) 

RUN(S,,,COMPILE, 

L4s0. 

7-8-9 (end-of-record, multiple-punched in column 1) 

*COPY, BMD3R 

7-8-9 (and-of-record, multiple-punched in column 1) 
(program control cards and data) 

6-7~8-9 (end-of-information, multiple-punched in 
column 1) 


For further details of the Source Library see the 
User's Manual, Section 5.1.2. SOURCE LIBRARY. 
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RAL PLOT INCLUDING HISTOGRAM 


GENc 


1. Program Name: BMD5D 
2. Central Memsy (CM}: 65000 
3. General Description 


a. This program provides a method by which graphs and histograms can 
be produced. 


b. Output for this program includes: 
1) GRAPHS. Two methods of plotting are available: 


a) This first method gives a one-page graph which has 50 units 
vertically and 100 units horizontally. The points are 
automatically scaled to conform te these dimensions, and 
3a scale is printed both horizontally and vertically. The 
points (data cards) need be in no special order. 


Say b) The second method gives a multiple-paae graph with as 
a many units vertically as there are values of the base 
variable. The values of the base variable (data cards) 
must be ordered and consecutive. The base variable is 
not scaled, The cross variables are scaled by the computer 
to conform to a horizontai dimension of 100 units. 


2) HISTOGRAMS 


A one-page histogram can be npraduced, with a maximum of 34 
intervals. The width of the interval must be specified; how 
ever, if the specifieu width would result in more than 34 
intervals, the program will print comments to this effect and eae 
will compute a new width which will give exactly 34 intervals. beac 
Scales are printed on the vertical and horizontal axes. 


c. Limitations per problem: 


1) p, number of original variables (1 < p< 500) 

2) on, number of cases (2 <n < 15000) 

3) 4, number of variables added to the original set after 
transgeneration (-499 < q < 499) 

4) p+q total number of variables (1 < ptq < 500) 

5) (p+q)n total number of data (2 < (prq)n < 15000 

6) Om, number of Transgeneration Cards (0 < m< 999) 

7) =k, number of Variable Format Cards (1 < k< 10) 


d. This program allows transgeneration. Codes 91-14 of the trans- 
generation list may be used. 


4. Order of Cards 


Cards indicated by letters enclosed in parentheses are optional. 
All other cards must be included in the order shown, 


| 


ef ce 
a. Job Card ## 2 IIIA 
(b.) REQUEST Card(s) *# - TIT B 
A c. LIBCOPY Card ee - ITI C 
d. LGO Card ## © III D) 
e. 7-8-9 Card (multiple-punched in column 1) 
f. Problem Card 
g. f-type Variable Format Card(s) ## - II B) 
(h.) Data Input Cards HH = ¥ B) | 
(Place data input deck here if data is from | 
cards) 
(i.} Standard Transgeneration Card({s) (## = II A) 
4 j. Selection Card ) 
k. Heading Card(s) te Repest from each 
| (1.}) Cross=Variable Card) | graph or histogram 
t 
{ Repeat f. through (1.) as desired 
; 
i . ee 
m. Finish Card (## - II D) 
| n, 6-7-8-9 Card ({multiple-punched in column 1) 
i 5. Preparation of Cards Specific for this Program 
i 
it f. Problem Card 


Col. 1-6 PROBLM 
7-12 Alphanumeric problem name 
13-15 Number of original variables (1 < p< 500) 
a 16-20 Number of cases (2 < n< 15000) 
. 21-23 Number of Selection Cards 
24-27. Number of variables added to original set after 
transgeneration (-499 < q < 499) 


AE Nee thereat an 


oe 


Note: (2 < (prq)n < 15000) 


LTTE egE MOR SE: MOT YR 


28-65 Leave blank. 
66-68 Number of Transgeneration Cards (0 < m< 999) 


5) To indicate for graphs how many variables are to be plotted 
against the base variable. (< 14) 

6) To indicate for graphs the choice of the type of graph. 

7) To indicate for histograms the width of an interval. 
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& 69-70 00 Data input from cards 
= 08 Data input from logical tape unit 8 
x 71-72 Number of Variable Format Cards (1 < k< 10) 
: ; j- Selection Card 
i A Selection Card has seven purposes: 
i i) Yo indicate whether a list of the data input is desired. : 
i 2) To indicate whether a graph or a histogram is to be produced. ee 
3) To indicate the base variable of the graph or histogram. # 
; 4) To indicate the number: of lines.of heading desired for each ate 
: i graph or histogram. a 
t 
t 
' 
J 
i 


If the Selection Card specifies that a graph is to be printed, the 
Heading Card is followed by a Cross-Variable Card which indicates 
the cross variables to be plotted against the base variable and the 
symbols used for each cross variable. 


Col. 1-6 = SELECT 
7 Number of lines in a heading. Each Heading Card 
specifies one line of printed output. The maximum 
number of lines allowed in the heading is two. (See 


_ card k.) 
5b 0 Listing of input data is not desired 
] Listing of input data is desired. 


9-10 Number of cross variables to appear on this graph : a 
{maximum is 14) 

tl-13. Index of the base variable. On graphs, the base 
variable will appear on the vertical axis. On i 
histograms, the base variable will appear on the 
horizontal axis. 

14-2 Form of the graph or width of interval if a histogram. 


Col. 1415 01 If a one-page graph is desired 
-1 If a multiple-page graph is desired, or 
Col. 14-24 Width of the interval for a histogram 
(punch the decimal point). 


k. Heading Card(s | 


Col. 1-72 Punch the desired heading. Each card is a line of 


the heading. There must be at lease one Heading Card, 


but no more than two, per graph or histogram. | 


(1.) Cross-Variable Card 


The Cross~Variable Card is punched as follows (for graphs only, 
not histograms), The cross variables specified to be crossed 
with one base variable will appear on one graph; the cross vari- 
ables will appear on the horizontal axis. 


Col. 1-6 CRSVAR 
7-9 Index of the Ist cross variable 
10 Symbol for the Ist cross variable (see below) 
11-15 Leave blank. 
16-18 Index of the 2nd. cross variable 
19 Symbol for the 2nd. cross variable 
20-24 Leave blank. | 


61-63 Index of the 7th cross variable « 
64 Symbol for the 7th cross variable 
65-69 = Leave blank. 


The symbols to be used for each cross variable must be specified. 
Allowable symbols are: 


»-,7 SJSKLMNOPQRSTUVWX YZ * 


The following symbols may not be used because they have been used to 
represent ties (more than one point occurring at the same coordinates) : 


Symbol No. of Points Symbo] No, cf Points Symbol No. of Points 


2 2 8 8 E 14 
3 3 9 9 F 15 
4 ky A 10 G 16 
5 5 8 VW H 17 
6 6 C 12 I 18 
7 7 D 3 / more than 18 


If there are more than seven cross variables, continue punching a 
second card in the same manner 


Col. 1-6 — CRSVAR 

7-9 Index for the 8th cross variable 
i 10 Symbol for the 8th cross variable 
1-15 Leave blank. 


61-63 Index for the 14th cross variable 
i €4 Symbol for the 14th cross variable 
65-69 Leave blank. 


The maximum number of cross variables for a specified base variable 
is 14, 


6. For 4 brief description of the computational procedure, refer to the 
Biomedical Computer Programs manual. 
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CORRELATION PROGRAM 


1. 
2 
3. 


4, 


Program Name: BMD2D 
Central Memory (CM): 75000 
General Description 


a. This program computes simple correlation coefficients, averages 
and measures of dispersion on entering variables and/or trans- 
generated variables. 


A special feature of this program is the selection of cases 
from the input data by specifying a boolean expression; i.e., 
case is accepted if it is in agreement with the expression; 
otherwise, the case is skipped. The expression consists of 
variables and constants involving relationships of equality or 
inequality written in a logical form using the operations 

AND and OR. 


b. Output from this program includes: 


1) Sums 

2) Means 

3) Standard deviations 
4) Correlation matrix 


Optional output includes: 


5) Cross-product deviations 

6) Variance-covariarnce matrix 

7) One-page cross-tabulation plots of any two variables, 
automatically scaled to 50 (vertical) by 100 (horizontal) 
character spaces or units. 


c. Limitations per problem: 


1) p, number of original variables (2 < p< 135) 

2) ©, number of original cases (2 < n< 99,999) 

3) Jj, number of Plot Selection Cards (O< < j < 99) 

4) q, number of variables added to the original set after 

transgeneration (-133 < q< 133) 

5) b, number of Case Selection Cards (O< b< 9) 

6) m, number of Transgeneration Cards (0 rem 
7) k, number of Variable Format Cards (I< k 


d. The program allows transgeneration of the input data. Codes 
01-16 and 41 of the transgeneration list may be used. 
Order of Card3. 
Cards indicated by letters enclosed in parentheses are optional. Al! 
other cards must be included in the order shown. 
a. Job Card a - IIL . 
(b.} REQUEST Card(s) +t - ITI B 
#4 = PREFACE to che CDC 6500 STATISTICAL PROGRAMS 
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e. 7-89 Card (multiple-punched ... column one) 
f. Problem Card 


c. LIBCOPY Card hes III 
J. L60 Cerd | (es TIT D | 
(g.) Standard Transgeneration Card(s) (## ~ IT A) | 
(h.) Case Selection Card(s) } 
i. Fetype Variable Format Card(s) ww = II B) 
(j.) Oats Input Card Me - 1B) 
(Place data input deck here if data 
is from cards) 
| (k.) Plot Selection Card(s) 
| Repeat f. through (k.) ss desired. 
i 
| 1. Finish Card (## - II D) 
i m. 67-89 Card (multiple-punched in column one) 
we 
5. Preparation of Cards Specific for this Program 
» | f. Problem Card 
Col. 1-6 — PROBLM 
\ 7-12 Alphanumeric problem name 
13-15 Number of original variables (2 < p< 135) 
16-20 Number of original cases (2 < n< 99999) 
21-22 Number of Plot Selection Cards; if none, 
i leave blank. (0 < j < 99) 
23-26 0000 No variables added to, or subtracced 
| from, the original set after trans- 
generation 
+q q variables added to the original set 
-_ CU after transgeneration (2 < p+q< 135) 
~ | -q q variables subtracted from the original 
ae set after transgeneration 
27-28 00 = No Case Sslection Cards 
+b b cards used for Boolean expression; 
case selection occurs after trans- 
generation (b < 9) 
“b b cards used for Boolean expression; 
— case selection occurs prior to trans- 
: generation (|b| < 9) 
29-39 NO if matrix of cross products is not | 
desired | 
31-32 NO if covariance matrix is not desired 
” 33~34 NO if aiternate input tape is not to 
be rewound 


35-65 Leave blank. 
66-68 000 No transgeneration 
m m Transgensration Cards (0 < m< 150) 
69-70 00 Data input from cards 
08 Data input from logical tape unit 8 
71-72 Number of Variable Format Cards (1 < k< 10) 
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(h.) Case Selection Card(s) 


ey Soares 


It is often useful to select cases if the value of a parti- 
cular variable is less than some constant, greater than 
some constant, equal to some constant, etc. Symbolically, 


where I is the index of some variable. To select only those 
cases where the values of a variable are between two con- 
stants involves the operation AND. 


v(I) > C AND V(I) < B 


{ 

| To select only those cases where either of two variables 
| must satisfy a relationship involves the operation OR. 
{ 


v(t) > C OR V(J) < B 
Perhaps a more complicated expression is desirable, e.g., 
(v(I) > A) OR (V(J)< B) AND (V(K) = Cc), ... 


According to rule. the entire Boolean expression is either 
true or false for the case being tested. It is examined 
from left to right. If an OR is encountered, and the 
| expression preceding the OR is true, the entire expression 
| is considered to be true for this case, and the case is se- 
lected for inclusion. 


Since parentheses cannot be used for compound AND/OR 
expressions, AND is assumed to precede OR. The statement 


. WOR X AND Y OR Z 
oN will operate as 
WOR (X AND Y) OR Z, 


A Case Selection Card is written as a sequence of con- 
ditions separated by an operation. A condition is a 
variable and a constant separated by a relationship. 


Variables: A variable is specified by the alphabetic V and 
the variable index; V(100), v(Ol0), v(149), v(008). The 
three-digit index is necessary; it is enclosed by perenthe- 
ses. 


F Constants. Constants are specified by their literal value, 
€.g., ~22.43, .99090, 1.0000, .00009. Five numeric 
characters with a decimal point are allowec. If the 
; sign (+,-) is used, then only four numeric characters 

7 are allowed. 
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Realetionships: Relationships are specified by using the 
following two-character codes: GT (grester .han), 
LT (less then), GE (greater than or equel to), LE 
than or equal to}, EQ (equal to}, RE (not equal 
tos. 


Gperations: Operations are specified by using the follow 
ing two-cheracter codes: AN (anc), OR (or), ** (end 
of expression), 


NOTE: (not greater than) -~—~ LE 
(not less than) -—- GE 
net greater than or equal to) =~ LT 
not less than or equal to) --—- GT 

Exemples: 


i) (V(OO2NEV (106) )** 


The case is accepted if variable 2 is not 
equal to variable 100. 


ii) . (V(230)GE100.00)AN(V(010)LT200,00) ** 
The case is accepted if variable J0 is greater 


than or equal te $00.00 and variable }0 is 
Tess than 200,00. 


The prepavation of the Case Selection Card is as 
foi lews: 


Col, 1-3 (Vi 
4-6 Three-digit variable index 
) 


&5 Two- character relationship 
10-11 v( 
12-14  Three-digit variable index 
{10+ 15 Constant (Keypunch decimal)} 
) 
17-18 Two~ character operstion 


This format is repcated four times per card endin 

in Column 72. The maximum number of cards is niné, 

The last operation of the expression must be **, There- 
fore, the user may specify from one to 36 conditions, 
each condition followed by an operation, the lest 
operation being **. 


(k.) Plot Selection Card(s) 


Col. 1-6 PLOTSL 
7-9 Index of the base variable 
(xX + axis) 
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69-71 


Number of variables to be cross-plotted with this 
base variable ( < 20) 

Index of the st variavle to be cross- 

plotted with this base variable 


Index of the and variable to be cress- 
plotted with this Lase variabie 


Index of the 20th variable to be cross- 
plotted with this base variable 


Each Plot Selection Card is independert. The same 
or d‘fferent base variables may be specified on 
additional cards, The maximum number of Plot Selection 


Cards is 99. 


6. For a brief description of the computational procedure, 
refer to the Biomedical Computer Programs Manua}, 
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